APR111919 
Oni. o¢ 


AMERICAN 


Journal of Mathematics 


EDITED BY 


FRANK MORLEY 


WITH THE COOPERATION OF 
A. COHEN, CHARLOTTE A. SCOTT, A. B. COBLE 


AND OTHER MATHEMATICIANS 


PUBLISHED UNDER THE AUSPICES OF THE JoHNS HopxKINs UNIVERSITY 


Ipaypdtwy Eleyyor ob Bhenopéevwr 


VOLUME XLI, NumsBer 2 


BALTIMORE: THE JOHNS HOPKINS PRESS 


LEMCKE & BUECHNER, New York. E. STEIGER & CO., New York. WILLIAM WESLEY & SON, London 
G. E. STECHERT & CO., New York. ARTHUR F. BIRD, London, A. HERMANN, Paris. 


1919 


Entered as Second-Class Matter at the Baltimore, Maryland, Postoffice. 
Acceptance for mailing at special rate of postage provided for in Section 1103, Act of October 3, 1917. 
Authorized on July 3, 1918. 


/ 


Asymptotic Satellites near the Straight-Line Equilibrium 
Points in the Problem of Three Bodies. 


By Danie, BucHANAN. 


§1. SaTELLITEs. 


It was shown by Lagrange in a prize memoir* in 1772 that if two finite 
spheres revolve in circles about their common centre of mass, then there are 
three points on the line joining their centres at which an infinitesimal body 
would remain if it were given an initial projection so as to be instantaneously 
fixed with respect to the moving bodies. These points are called the straight- 
line equilibrium points of the problem of three bodies. Starting from minus 
infinity the order of the equilibrium points and the finite bodies is an 
equilibrium point, a finite body, a second equilibrium point, the other finite 
body, the third equilibrium point. 

If the infinitesimal body is given an initial displacement from a point of 
equilibrium and initial conditions are so chosen that it moves in a closed orbit 
relatively to the moving system, it is then called an oscillating satellite. 

The problem of the oscillating satellite has been discussed extensively. In 
the papers cited below,+ the differential equations are limited to their linear 
terms and the orbits restricted to the plane of motion of the finite bodies. 

A rigorous demonstration for the existence of periodic orbits for the 
oscillating satellite and a practical method for constructing them are given by 
Moulton in Chapter V of his “Periodic Orbits.”t In this memoir the differ- 
ential equations are unrestricted in the number of terms which may be taken, 


* Lagrange, “Collected Works,” Vol. VI, pp. 229-324. 

+ The following references to the literature of the oscillating satellite are taken from Chapter V 
of Moulton’s “ Periodic Orbits”: 

Poincaré, Les Méthodes Nouvelles de la Mécanique Céleste, Vol. I (1892), p. 159. 

Burrau, Astronomische Nachrichten, Nos. 3230, 3251 (1894). 

Perchot and Mascart, Bulletin Astronomique, Vol. XII (1895), p. 329. Moulton adds, “apparently 
their work is vitiated by an error in establishing the existence of the solutions, and their construction 
fails where they stopped.” 

Sir George H. Darwin, Acta Mathematica, Vol. XXI (1897), p. 99. 

Plummer, Monthly Notices, Royal Astronomical Society, Vol. LXIII (1903), p. 436, and Vol. LXIV 
(1903), p. 98. 

{This memoir will be cited as the “Oscillating Satellite.” Another method is given in Chapter 
VI of the “ Periodic Orbits,” but we are concerned with the former method only. 
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and the orbits are not limited to two dimensions as in the previous literature. 
Three classes of orbits are shown to exist and they are designated according to 
their periods as Class A, Class B, and Class C. The orbits of Class A and 
Class C are of three dimensions, while those of Class B are of two dimensions. 
The orbits of Class C are shown to exist under special conditions which might 
never be realized in the problem. Their period is a multiple of the periods of 
Class A and Class B when these latter periods are commensurable. Practical 
constructions are made for the orbits of Class A and of Class B, but, owing 
to the complexity of the problem, no attempt has been made to determine 
whether orbits of Class C exist which are distinct from those of Class 


A and Class B. 


§2. Asymptotic SaTELLITEs. 


The object of this paper is to obtain solutions of the differential equations 
of motion of the infinitesimal body which will approach the periodic solutions 
of Class A and Class B as the time approaches plus infinity or minus 
infinity. When the infinitesimal body moves in an orbit defined by such 
solutions, it will be called an asymptotic satellite. 

The question of the existence of solutions which are asymptotic to the 
orbits of Class C is not considered in this paper, owing to the fact that it has 
not been determined in the “Oscillating Satellite” that orbits of Class C exist 
which are distinct from those of Class A and Class B. 

The orbits which are asymptotic to the equilibrium points themselves have 
been determined by Warren.* These orbits are of two dimensions and lie in 
the plane of motion of the finite bodies. 

The form of the asymptotic solutions is that adopted by Poincaré,t each 
term being of the type e’P(t), where 4 is a constant having its real part 
different from zero, and P(t) is a constant or periodic function of ¢. It has 
been shown by Poincaré that solutions of this type will converge for all values 
of ¢, provided that certain divisors which appear in the construction of such 
solutions are different from zero.{ If these divisors vanish, terms of the form 
te“P(t) will arise. If, therefore, the construction can be made so that no 
terms occur in ¢ explicitly, the divisors previously mentioned are different from 
zero and the solutions will converge for all values of ¢. Hence it is sufficient 


* Warren, “ A Class of Asymptotic Orbits in the Problem of Three Bodies,” AMERICAN JOUBNAL OF 
MatTuHeEmartics, Vol. XXXVIII, No. 3 (1916), pp. 221-248. 

¢ Poincaré, “ Mécanique Celéste,” Vol. I, p. 340. 

¢ Poincaré, loc. cit. p. 341. 
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to consider only the formal construction of the asymptotic solutions, and if 
solutions can be constructed so as to contain no terms in ¢ explicitly, their 
convergence is assured by Poincaré’s theorem. 


§3. Tue DirrerentiaL Equations or Motion. 


Let the motion of the infinitesimal body be referred to a set of rotating 
rectangular coordinates £y¢ of which the origin is at the centre of mass of the 
finite bodies, the £-axis is the line joining the finite bodies, and the &y-plane is 
the plane of their motion. The &- and y-axes rotate about the ¢-axis in the 
direction of the motion of the finite bodies and with the same angular velocity. 
The units of length, mass, and time will be taken so that the distance between 
the finite bodies, the sum of their masses, and the Gaussian constant respectively 
shall each be unity. With the units thus chosen, the mean angular motion of 
the system is likewise unity. Let the masses of the finite bodies be denoted by 
1—y and uw, 0<u<4. On denoting the coordinates of the infinitesimal body 
by &, , ¢, and differentiation with respect to ¢ by primes, the differential 
equations of motion are* 


dU dU oU ) 


2(1— 2 


= (14) (18+ =) 
n= +h. | 
The points of equilibrium are the solutions of the equations t 


0. 


There are two sets of points which satisfy these equations. One set consists 
of the two points which are at the vertices of the two equilateral triangles on 
the opposite sides of the line joining the finite bodies. The orbits which are 
asymptotic to these points and also to the periodic oscillations near these 
points are discussed in another paper.{ The other set of points consists of 


* Moulton, “Introduction to Celestial Mechanics” (1914), p. 279. 

+ Moulton, “ Introduction to Celestial Mechanics,” p. 290; Charlier, “ Die Mechanik des Himmels,”’ 
Vol. II, pp. 102-111. 

{ This paper is now under consideration for publication. 
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three points which lie on the straight line joining the finite bodies. Let the 
coordinates of these points be denoted by &, 0,0, where the particular value 
of & depends upon the equilibrium point in question.* The points themselves 
will be denoted by (a), (b), and (c), where (a) lies between +o and the 


finite mass u, (b) between uw and 1—y, and (c) between 1—yu and —o. 
If the infinitesimal] body is given a small displacement from an equilibrium 


point and a small velocity with respect to the finite masses such that 
¢=0+2, (2) 
+2’, n'=0+Y’, ¢=0 +2, 

then the differential equations (1) become t 


—2y’ = = +X 2), 


dU 
+22 (Z, 2) 
Oz =22 (%, 2), 
2 2 ‘i 
U =4(1—u) >) 
1 2 
where X, Y, and Z are power series in %, y’, 2. These series converge within 
certain regions about the equilibrium points. + 


§4. THe Pertopic Orsits. 


In showing the existence of periodic solutions of the differential equations 
of the “ Oscillating Satellite” which correspond to equations (3), and later in 
making the construction of these solutions, the transformations 

Yoey, Z=ez, (4) 
are made, where ¢ is an arbitrary parameter and 6 is determined as a function 
of e so that the solutions in x, y, and @ shall be periodic with the periodic 
2x in t. On denoting differentiation with respect to t by a dot over the 
variables, the differential equations (3) become as a consequence of (4)§ 

e+....+X,e*1+....], 
(148) Yet... (5) 
| 


*< Oscillating Satellite,” equations (4). ¢ “Oscillating Satellite,” equations (6). 
¢ “ Oscillating Satellite,” §77. § “ Oscillating Satellite,” equations (11). 
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where X,, Y,, and Z, are homogeneous polynomials of degree k in 2, y, and z. 
From (3) it is obvious that the X, are even in y and z, Y, odd in y and even 
in 2, and Z, even in y and odd in z. The explicit values of these terms up to 
k=3 are* 
=(1+2A)a, X,;=2C 
Y,=(1—A)y, Y.=3Bay, Y3;=3C 
Z,=—Az2, 2,=3Bzz, Z,=4C (—4e2+ y'%2+2), 


1— 


+1—n+ 


C= 


The upper, middle, or lower signs are to be taken in B according as the 
equilibrium point is (a), (b), or (c) respectively. 

If periodic solutions of (5) exist, their periods are determined from the 
periods of the solutions of the linear terms of (5). The solutions of these 


linear terms are t+ 


a= i=V—1, 

2=K,cos VAt+K, sin V At, 


where K,,...., K, are the constants of integration, and o* and p’ are the 
negative and positive roots respectively of the quadratic in 4’, 

A‘ + (2—A)A?+ (1—A) (14+2A4) =0. (7) 
There are three real periods in these solutions, viz., 22/ VA, 2n/c, and 
P=2jn/VA=2kn/c, where j and k are positive integers and VA and c are 
commensurable.t{ These are the periods of the orbits of Class A, Class B, 


and Class C respectively. 
Orbits of Class A. 


The periodic solutions of Class A are§ 
= =0e+ (a,+b, cos 2VAt)2+0E+...., 
= ey, =0e+ (c, sin 2V At)? +0F+...., 


2, sin Var) e+0e+d,(3 sin VAt— sin3 VAt)eé+...., (8) 


§ =8,=0e+ bP P+... J 


* «< Oscillating Satellite,” equations (53), (15), (35). 

7 “Oscillating Satellite,” equations (16), (27), and (28). 

tIn §83 of the “Oscillating Satellite” it is shown that there are infinitely many values of x 
between 0 and 4 such that VA and o are commensurable. 

§ “ Oscillating Satellite,” equations (71). 
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where —3B 3B(1+3A) 
4A4(1+2A)’ 44 (1—7A +184?) ’ 
644°*| 1—74+ 184° 
9 3B*(1—34 +144?) 
16A? L (1+ 2A) (1—74+ 18?) 


The expressions for %, and Y, are power series in ¢’ with sums of cosines 
and sines respectively of even multiples of V At in the coefficients. The highest 
multiple of V.At which appears in the coefficient of e* in each seriesis k. The 
solution for Z, is a power series in odd powers of ¢ having sums of sines of 
odd multiples of V.At in the coefficients. The highest multiple of VAt which 
occurs in the coefficient of e* is likewise k. The initia] conditions are so chosen 
in the construction of these solutions that 


z(0)=0, 2(0) =e, 


that is, the parameter « is proportional to the initial projection of the infini- 

tesimal body from the plane of motion. The actual projection is e/(1+4,). 
Numerical examples of these orbits have been considered in the “ Oscillating 

Satellite” and the following results have been obtained when the ratio of the 


finite masses is ten to one, or 1—w=10/11, w=1/11." 


TABLE 1. 

Coefficient. Point (a). Point (b). Point (c). 
A 2.548 6.510 1.082 
o? 2.811 6.820 1.144 
p” 3.359 11.330 0.226 
n 2.657 3.990 2.014 
m —0.747 — 0.397 —3.091 
B 6.548 —10.961 —1.136 
C 18.283 090.740 1.196 
a; —0.316 0.090 0.249 
b, 0.151 | — 0.086 —0.230 
—0.112 0.018 0.226 
d, 6.00 —0.002 
6? 0.184 0.467 0.001 


* This was the ratio used by Darwin in the memoir cited at the beginning of this article. The same 
ratio was adopted by Moulton in the “Oscillating Satellites” so that he might compare his results with 


those of Darwin. 
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Orbits of Class B. 
The periodic solutions of Class B are * 

eX, = (cos ot)e+ (a,+b, cos ot +0, cos 2at)e?+...., 

(—nsinot)e+ (—nb,sin ot +d, sin 2ct)e’+...., (9) 

6 
where a, b,, c,, ...., are known constants. The algebraic expressions 
for these constants may be found in equations (101), (103), and (107) of the 
“Oscillating Satellite.” Their particular values for 1—u=10/11, uw=1/11 are 
to be found in Table 2. The first seven coefficients which appear in Table 1 
have been used in determining the values in Table 2. 


TABLE 2. 
Coefficient. Point (a). Point (0). Point (c). 
le —4.078 8.172 0.554 
b, 1.996 — 4.976 —0.038 
Ce 2.082 — 3.196 —0.516 
nb, 5.305 | —19.855 —0.077 
d, 1.250 — 1.478 —0.276 
6? 3.955 8.553 —1.407 


The expressions for % and ¥, in (9) are power series in e with sums of 
cosines and sines respectively of ot in the coefficients. The highest multiple 
of ot which occurs in the coefficient of e* is k. The initial conditions are so 
chosen for the construction of these solutions that 


%(0)=«, Yy.(0)=0. 
In this case the parameter ¢ denotes the initial displacement along the £-axis 
of the infinitesimal body from an equilibrium point. 


§5. Tue Equations or VARIATION. 


Suppose now that the initial conditions for the periodic orbits are 
changed slightly. Then the infinitesimal body will deviate from its periodic 
orbit, and the amount and character of the deviation will depend upon the 
initial conditions chosen. The equations for the disturbed orbit, in terms of 
the variables x, y, and 2, will be 


y=yjtq, 2=24+1;1; (10) 


* “ Oscillating Satellite,” equations (97), (105). 
+The r, and r, in this equation and in subsequent equations are different from r, and r, in 
equations (1) and (3). 
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where ex;, ey;, and ez; represent the periodic orbits of Class A or Class B 
according as j=1 or 2 respectively. The additional terms p,, g;, and 1; are 
functions of t which depend upon the initial disturbance from the periodic 
orbit. We desire to determine these functions so that they shall have the 


asymptotic form already referred to in § 2. 
On substituting (10) in (5), the differential equations which define p,, q;, 


and r; are found to be 


(1+6,)°P;, 
Gj+2(1+4;) (1+4;)? = (1+8,)'Q;, (11) 
(1 +R = (1+6;)*R,, 


where 
P,=—1—2A+6Bex,—6Ce’ 


e=A—3Bex,+3Ce +...., 
P; 
Q; 
As we propose later to construct numerical solutions which are asymptotic 
to the particular periodic orbits mentioned above, we insert here the particular 


values of the various P,;,, Q,,, and R,, which are obtained by making use of the 
numerical values listed in Tables 1 and.2. The following results are found. 


Orbits of Class A. 


Equilibrium point (a). 
Py = —6.096 +2? (8.384—15.686 cos 2V.At) + ()et+...., 
Qa =e? (2.200 sin 2V At) + ()et+...., 
Py3= Ry = —e (12.297 sin VAt) + ()e+...., 
cos 2VAt)+()et+...., 
Ry=2.548—e(9.918—13.160 cos 2V At) +()e'+..... 


Equilibrium Points in the Problem of Three Bodies. 


Equilibrium point (b). 
.328 cos 2V At) +()ett+...., 
(0.592 sin2VAt) + ()ett+...., 
=R,, =e (12.890 sin VAt) + ()P+...., 
Ry cos 2V At) +()ett+..... 
Equilibrium point (c). 
Py=—3.164+2(1.616—1.746 
Pye=Q =e (0.770 sin 2V-At) + ()et+...., 
Py=Ry=e(3.275 sin VAt) + ()e+...., 
cos 2VAt) + ()ef+...., 
cos2VAt)+..... 


Orbits of Class B. 

Equilibrium point (a). 

cosot)+...., Po=Qy=e(52.27 sinot)+...., 
Qoo=1.55—e(19 . 64 cos oT) 
R,,=2.55—e(19.64 cosot)+..... 

Equilibrium point (b). 

Py=—14.02—e(65.77 coset) +...., 
cosot)+...., 
R. »=6.51+e(32.88cos ot) +..... 

Equilibrium point (c). 
Py=—3.16—e(6.82 cosot)+...., Poe=Qu=e(6.86 sinot)+...., 
cos ot) +. 
cosot)+..... 

If we consider only the terms of (11) which are linear in p,, q;, and 1,, 
we obtain the equations of variation.* For the orbits of Class A they are 


P.—2(1+4,) =0, 
Q+2(1+6,) (1481)? + =0, (12A) 
Bsn) =0, 
and for orbits of Class B, they are 
Qo +2(1+6,) Pot Qe292) =0, (12B) 


* Poincaré, loc. cit., Vol. I, Chap. 4. 


| | 
87 
| ‘ 
12 
{ 


88 Bucuanan: Asymptotic Satellites near the Straight-Line 


The periodic solutions (8) and (9) are called the generating solutions* of 
(12A) and (12B) respectively. 


§6. Sotutions or THE Equations or Variation (12A). 


The equations of variation are linear differential equations having periodic 
coefficients, the periods being the same as those of the corresponding generating 
solutions. Differential equations of this type were first treated by Hillt in 
his celebrated memoir on the lunar theory. A very extensive list of references 
to the literature of these differential equations was given by Baker in 
a recent memoir. { | 

The method which we shall adopt in the construction of the solutions of 
these equations is the one developed by Moulton and Macmillan.§ The form]| 
of the solution is, in general 

(7), 
where (7) is a periodic function having the same period as the coefficients 
of the differential equations, and a is a constant called the character- 
istic exponent. 

It has been shown by Poincaré that if the generating solutions contain 
an arbitrary constant which does not appear explicitly in the original 
differential equations, then a set of solutions of the equations of variation is 
obtained by taking the first partial derivatives of the generating solutions with 
respect to this constant. The characteristic exponent associated with this set 
of solutions is zero, and the solutions themselves are either periodic or consist 
of periodic functions plus t+ times other periodic functions. The initial 
time ¢, is one such arbitrary constant which is present in all differential 
equations of dynamics in which the potential function e.g., U in equations (1), 
does not contain the time explicitly. Differentiation of the generating solu- 
tions with respect to this constant yields a set of periodic solutions. Another 
such constant usually present is the scale factor of the generating solutions. 
Differentiation of the generating solutions with respect to this constant yields 
a set of solutions comprised of a periodic function plus Kt, K a constant, 
times the periodic function obtained by differentiating with respect to ¢,. 


* Poincaré, loc. cit. Vol. I, Chap. 4. 

+ G. W. Hill, “ Collected Works,” Vol. I, p. 243; Acta Mathematica, Vol. VIII, pp. 1-36. 

¢H. F. Baker, “On Certain Linear Differential Equations of Astronomical Interest,” Philosophical 
Transactions of the Royal Society of London, Series A (1916), Vol. 216, pp. 129-186. 

§ Moulton and Macmillan, “On the Solutions of Certain Types of Linear Differential Equations 
with Periodic Coefficients,” AMERICAN JOURNAL oF MATHEMATICS, Vol. XXXIII (1911), No. 1, pp. 63-97. 

|| Floquet, Annales Scientifiques de V Ecole Superieure, Series II, Vol. XII (1883), p- 47. 

{{ Poincaré, loc cit,, Vol. I Chap. 4. 
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Let us consider first equations (12A), and those solutions which have 
characteristic exponents different from zero. These equations are simul- 
taneous and must be considered together. To find their solutions let 

(13) 
where K,, Z,, and M, are arbitrary constants. An existence proof, which is 
omitted here, would show that 


a, =a taM 

=U (14) 
= VO M P+... 


where the af" are constants so determined that the various u?, v@, and 
w+) shall be periodic with the period 2%/ VA int. On substituting (13) in 
(12A) and considering only the terms independent of e¢, it is found that af 
must satisfy the biquadratic (7) if the arbitrary constants K,, L,, and M, are 
to have values other than the trivial ones K,=L,—M,=0. It is shown in the 
“ Qscillating Satellite,” § 82, that this biquadratic admits two purely imaginary 
solutions, denoted by +o7 and —ci, and two real solutions, denoted by +p 
and —p, for all values of uw such that 0<yu<#, and for each equilibrium point 
(a), (b), and (c). 

Since the functions u,, v,, and w, are multiplied by arbitrary constants, 

we may assume without loss of generality, that u,(0)=1. Hence 

u(0)=1, u?”(0)=0 (Ba: . (15) 
Then v,(0) and w,(0) can be determined from the differential equations which 
define v, and w,. 

The various solutions of (12A) may be constructed by substituting (14) 
and (13) in (12A), and giving to a{” the values oi, —ci, p, and —p in turn. 
By equating the coefficients of the various powers of ¢ we obtain sets of differ- 
ential equations which can be integrated step by step. The arbitrary constants 
of integration arising at each step can be uniquely determined from the condi- 
tions (15). The various af” can likewise be uniquely determined by the 
condition that v@, and shall be periodic with the period 
int. From Poincaré’s extension to Cauchy’s theorem,* it is known that these 
solutions will converge for ¢ sufficiently small numerically and for all values 
of t such that 0<7t<Z, where 7 is an arbitrary period chosen in advance. 
Obviously, we may choose the period 22/ VA. 


* Poincaré, loc. cit., Vol. I, p. 58. 
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Let us denote the solutions of (12A) when af=«ai by 

Uy, wy. (16) 
The arbitrary constants are omitted here but they appear in equations (25). 
It is found that o, is a power series in e with real constant coefficients. The 
functions u,, and v,, are likewise power series in <’, but their coefficients are 
sums of cosines and i times sines of even multiples of VAt. The function w,, 
is a power series in odd powers of e, vanishing with e, and the coefficients are 
sums of cosines and i times sines of odd multiples of VAt. In each of these 
three functions, the highest multiple of VAt which appears in the coefficient 
of e* is k. So far as the computation has been carried out, we have 


Hy =1+()P+...., tu=nt()e+...., 


= 08 Var— Vaz |e+ Oe+...., 
—(8 + 7A—224?)' 8(1—A) 


Since the differential equations (12A) do not contain 1, they are unaltered 
by a change in the sign of 7. If then we change the sign of 1 in the differential 
equations and in the initial conditions (15), we obtain solutions which differ 
from (16) only in the sign of 7. Let this set of solutions be denoted by 


"Uy, Wy, (17) 
where U., Vz, and Wy» differ from uy, vy, and wy, respectively only in the 
sign of i. The same result would be obtained by putting a =—ie and pro- 


ceeding as in the construction of the solution (16). 
Likewise by putting a{ equal, first to p, and then to —p, we obtain the 


respective solutions 


and 


where p, is a power series in e” with real constant coefficients, Further us, uy; 
Vis, Vu; Wis, Wy are similar to uw, Vy, and w, respectively, except that the 
coefficients of the former functions are all real; also 
Uyg(T) =Uyg(—T), Wy (tT) 
The computation for the first terms of these series gives 
Ug 


3B 2 A: | 8 
cos VAt 


ly 
+3+7A—22A?} 4 6 (1—A) (1 |+ 
3 


9 
2 
| (44 +6") 
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Obviously, the solutions (19) differ from (18) only in the sign of 7, q,, 
and r,. This property can be derived directly from the differential equations, 


which we proceed to show. We shall show that if 
a(t), G=Fa(t), =F (20) 
are solutions of the differential equations (11), then 
Dj=Fy(—t), GQ=—Pp(—t), 
are likewise solutions. If this property holds for equations (11) it will 
evidently hold for equations (12A) or (12B), i. e., when the right members of 
(11) are omitted. 


Let us consider the differential equations (1). Since U is even in y and @, 
it follows that oe is even in and is odd in y and even in and 
is even in 7 and odd in @. After the substitutions (2) are made, the right 
members of (3) are functions of ¥ and Z as indicated. As the substitutions (4) 
do not alter the parity of y and 2, then it follows that the right members of 
(5) have the same parity in y and ¢z that the corresponding right members 
have in ¥ and Z. When the substitutions (10) are made in (5) we obtain 
(11), and the terms of (11) which carry (1+)? as a factor in each equation 
arise from the corresponding right members of (5). Let these terms be 
denoted by (11,1), (11,2), and (11,3) respectively. Then they have the 
following properties: 

(11, 1) is even in and even in {2,, 75}, 

(11, 2) is odd in {y,, q;}, and even in jz,, 1;{, 

(11, 3) is even in {y,, q;}, and odd in }z,, 7;}, 
where the braces {| denote that the variables within are to be considered 
together. Further, since 7 is a factor of the right member of the second 
equation in (3), it follows that (11,2) carries a factor which is odd in {y,, q,}. 
Similarly, since 2 is a factor of the right members of the third equation in 
(3), then (11, 3) carries a factor which is odd in }z2,, q;}. 

Now, if the signs of y,, 2;, g;, and r; be changed in the above expressions, 
then (11,1) remains unchanged while (11, 2) and (11,3) change signs. On 
examining the generating solutions (8) and (9) we observe that changing the 
signs of y; and 2; is equivalent to changing the sign of t in these solutions. 
Thus, if we change the signs of t, g;, and r;, (11, 1) remains unchanged while 
(11, 2) and (11, 3) change signs. If these changes are made in all the terms 
of (11), the first equation remains unchanged while each term of the remaining 
equations changes sign and the factor —1 can be divided out. Thus the 
equations (11) remain unchanged if we change the signs of t, g, and 1;. 


| 
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Suppose now that (20) is a set of solutions of (11). Let 


G=—-Q, T=—T. (21) 


Then the differential equations obtained by making these substitutions in (11) 
are identically the same in p,, g,, 7;, and 7 as (11) are in the former variables. 
With the same initial conditions for p,, g;, 7;, and tT as for equations (11), we 
obtain the solutions 

p;=F (7), G;=F p(T), 7;=F (7), 
where F',, F,, and F', are the same functions as in (20). On restoring the 
former variables by (21) we have 


g=—Fa(—t), 

as solutions of (11). Therefore if a set of solutions of (11) is known, another 
set can be obtained by changing the signs of +, g;, and r; in the former set. 

Let us next consider the two remaining sets of solutions of (12A). These 
solutions have characteristic exponents zero and may be obtained, as already 
referred to in this section, by differentiating the generating solutions (8) with 
respect to the arbitrary constants t, and e. Taking the constant ¢, first, we 
have the solutions 


0%, O02, Or 1 


1 


02, 0% Ot 


1 
n= =— | (cos VAt)yet | 


As these solutions are later multiplied by arbitrary constants, we may neglect 
1 
the factor — T+h and consider 
Pi=%55 H=%5, 1=Wis (22) 
as the solutions. Similarly, on differentiating (8) with respect to e we obtain 


_t | _ 


Since t= (t—t,)/(1+4,) and 6, is a function of ¢, then the constant ¢ enters 


| 
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the generating solutions not only explicitly but implicitly through 6, and 7. 
Hence 
\ de Ot 00, de’ Oe .\ de Ot 00, de’ 
_ (th) , Bh 
Oe \ de Ot 00, Oe’ 
where the parentheses around the partial derivatives denote that the differ- 
entiation is performed only in so far as « occurs explicitly. When the 
differentiations are performed, we have 


where 
= cos 2VAt)et+()P+...., 
ou 
Vy = ($2) (6, sin 
oz 1 


1 


It remains now to show that the solutions which we have obtained, viz. | 
(16), (17), (18), (19), (22), and (23) constitute a fundamental set. The 
criterion for a fundamental set is that the determinant formed from these 
solutions and their first derivatives with respect to t shall be different from 
zero. This determinant is 


(24) 


(—iw,s), Wis (—wy), W153, Wye t 


It is a constant,* and its value can be determined with the least difficulty by 
putting t=0. Thus we obtain 

A,=4ie(mp+nc) (mo—no) + terms of higher degree in e. 
It is shown in the “Oscillating Satellite,” equations (36), that mp+mno and 
mo—np are different from zero. Hence the determinant A, is different from 


* Moulton, “ Periodic Orbits,” Chap. 1, Sec. 18. 


| | 
| 
| | 
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zero for e not zero, but sufficiently small numerically. Therefore the solutions 
which have been obtained for (12A) constitute a fundamental set, and the 
most general solutions are 


Kye", + + Kye” + Kye Uy, + K 
+K6(uy 

Kye"vy, ivy, + —K vy, + K 
+K705), 

1, = iwy— Ke iw, + Kye" —K Wy 


where K,,...., K, are arbitrary constants. 


§7. THe Sotutions oF THE Equations or Variation (12B). 


The solutions of the equations (12B) are obtained in the same way as the 
preceding solutions were found. The construction is simplified by the fact 
that the last equation is independent of the first two. 

The two sets of solutions of the first two equations of (12B) which have 
characteristic exponents different from zero, are 


(26) 


and 


where p, is a power series in e with real coefficients, the terms of lower degrees 
being 

The functions w2,, V2, Weg, ANG Veg are power series in e with sums of sines and 
cosines of multiples of ct in the coefficients, the highest multiple of ot which 
occurs in the coefficient of e* being k. All the numerical coefficients in these 
series are real. From the property of the differential equations (11) which 
was proved in the preceding section, it follows that the solutions (27) are 
obtained from (26) by changing the signs of t and q,. Thus uw. and vz differ 
from u,; and v., respectively only in the signs of the sines. 

Since these solutions (26) and (27) are multiplied later by arbitrary 
constants, we may take u,.,(0)=1. The algebraic forms of the functions w,, 
and v,; are 


1 + (aq, cos ot Sin ot)e+()e’+...., 
= M+ (doy COS OT-+ Sin 
(T) (—T), =Vq,(—T). 


| 
| | 
| 

| 
| 
i 
| 
| 


Equilibrium Points in the Problem of Three Bodies. 95 


The two remaining solutions of the first two equations of (12B) are 
obtained by differentiating the generating solutions (9) with respect to ¢, and e. 
Differentiation with respect to ¢, gives 


_% % ~_% 1 &, 


and, since the multiplier — — may be dropped, the solutions become 
2 


= sin ot) e— sin ot + sin 2ot)e?+...., 
(28) 
Vy = = cos ot) e— (nb,o cos ct — cos 2Zot)e’+..... 
Differentiation with respect to « gives 
t+ Log, (29) 


where 
Uy = = cos ot +2 (a,+b, cos ct +, cos 2ot)e+.. 


sin ot —2 (nb, sin ot—d, sin 2oT)e+...., 
1 
1+6, 

The determinant of these solutions (26), (27), (28), and (29) and their 
first derivatives with respect to 7 is 


L=— 


(— Vee), Veg, Vat 
%), (P2V22—V22) ’ +L (TV 28), 


= 2ce(mp+nc) (mo—np) + terms of higher degree in e. 


As in the former determinant A,, equation (24), the factors mp+no and 
ma—np are different from zero, and since «0 it follows that A, +0 for e not 
zero, but sufficiently small numerically. Hence the solutions (26), (27), (28), 
and (29) constitute a fundamental set of solutions of the first two equations, 


and the most general solutions are 


= + Ugg + Lg + Ly (tog + , (31) 
do = Lye"? + 03 + Lg (Vou 
where L,,...., L, are arbitrary constants. 


13 
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The general solution of the last equation of (12B) is readily obtained, 
and it has the form 

= + Lee (32) 

where L; and [, are arbitrary constants, and ow, w.,, and We are power series 


in ¢ of the form 
3Be 


Wy=1+ cos ot—2iV A +...., 
3B 
Wee = 1+ dA) (o? cos ot sin ot—o?) +...., 


o= VA+0e+are?+..... 


Since LZ; and LZ, are arbitrary, the initial values of w.,; and w., may be chosen 
so that w.,(0) =w..(0)=1. The coefficients of the various powers of ¢ in these 
functions are sums of cosines and 7 times sines of multiples of ot, the highest 
multiple which occurs in the coefficient of «* being k. Further, w,, and wy, 
differ only in the sign of 7. In the series for a, all the coefficients are real. 
The determinant of the two solutions in (32) and their derivatives with 


respect to 7 is 


+ Wy), (—toWy + We), 
= —2ia+ terms in 


As 


and this is different from zero for |e| sufficiently small. Hence the solutions 
(32) constitute a fundamental set of solutions of the last equation of (12B). 


§8. THe Construction or Asymptotic SoLuTIONs. 


Having determined the solutions of the equations of variation, we shall 
now proceed to construct the asymptotic solutions of (11). We shall consider 
first the solutions which approach zero as t approaches +o, and then show 
how to obtain from these the solutions which approach zero as +t approaches 
— 2, 

In making the constructions, it is convenient to introduce a parameter y 
by the substitutions 


Di=DiY, (33) 
where p;, g;, and 7; are new dependent variables. Then, since the asymptotic 
solutions are to be constructed as power series in y, we put 


TOP Y TOP YH tapytaPy’ t+ (34) 


| | 
— 

| 
| 
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When (34) and (33) are substituted in (11), the factor y will divide out. Let 
the resulting differential equations be denoted by (11’). On equating the 
coefficients of the various powers of y in these equations, we obtain sets of 
differential equations which define the functions p{”, g®, and r§” in (34). In 
order that the solutions shall be asymptotic we must impose the condition 
(C,), that each term of the solutions for the various p™, g{, and r shall 
contain e~*” as a factor. 
As an arbitrary constant arises at each step of the integration we may impose 
the further condition 


(C,), that p;(0) 
from which it follows that 


(0)=a,, (k=1,....0). 


Asymptotic Orbits of Class A. 


Let us consider first the orbits which are asymptotic to the periodic orbits 
of Class A. We put j7=1 and consider the various differential equations 
obtained by equating the coefficients of the same powers of y in (11’). 

For the terms in p{, g&, and r® we obtain a set of differential equations 
which are the same as (12A), except for the superscript 0. The general solu- 
tions are the same as (25), and when the condition (C,) is imposed all the 
constants of integration must be zero except the one associated with e~”. 
From condition (C,) it follows that this constant must have the value a,. 
Therefore the desired solutions are 


pO=ae"u,, rO=—aje "wy. (35) 


The differential equations which are obtained from equating the coefficients 
of y to the first degree in (11’) have the same left members as (12A), except 
for the superscript 1 on the variables. On denoting the corresponding right 
members by P®, Q{, and R{®, we find that 


where U®, V&, and WY are power series of the same form as uy, Vy, and Wy 


respectively. The complementary functions of these differential equations 
are the same as (25), but let us suppose that the arbitrary multipliers are 


| 
| | 
| | 
| | 
| | 
| | 
| 
| | 
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kf, ....,k@. The particular integrals can be obtained by varying these 
multipliers, Thus 
=0, 
ke" + Pe" + (Pitas + 4s) 
+hP Krv45) =0, 
+ ( Vu) + + k@ [Viet K ] = 
=0, 
ke ef" (—9,Wy+iwy) + ke + ken" ( Ws) 
+ (P1Wu—Wy) + ke [Wet K (Wy +7 ] =RY. 


| (36) 


The determinant of the coefficients of k,...., kK is A,, equation (24), and 
since it is different from zero for e not zero but sufficiently small numerically 


these equations can be solved for k®,...., kK. Thus 
kP = (fuel, ...+¢@), (37) 


where A{? is the determinant formed by replacing the elements of the J-th 
column of A,, with 0, P®, 0, 0, and R{& respectively. Since and 
R contain no terms in e+” or e+", the integrations of (37) for kf, ...., 42 
will yield no terms in ¢ explicitly. Such terms will occur, however, in the 
integrations for kf and kf, but when they are substituted in the complementary 
functions they will cancel off. The complete solutions are thus found to be 


=K Pe Pe + Wy —K Pe Wy 
1 2,72, 1 
+K§? (Wy 


r (38) 


J 
(1) 


where K{,...., K@ are the constants of integration, and and 
are functions similar to uy, vy, and wy respectively. On imposing condition 
(C,) upon these solutions we obtain 


| 
} 
a 
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and from condition (C,) we have the remaining constant 

KP (0). 
When the constants of integration are thus determined, the solutions (38) 
take the form 


=ea2(e +e" wD), 


where the functions u{?, uf; v{P, vo; wh, wY are of the same form as u,, 
Vy, and respectively. 

The remaining steps of the integration can be carried on in essentially 
the same way. By an induction to the general term we shall show that the 
solutions of (11) can be constructed to any desired degree of accuracy. 

Let us assume that p{”, qi”, r@ have been constructed for y=1,2,....,4—1, 
and that 

y+1 
where the functions u{?, v{P, and w{? are similar to u,,, vy, and wy, respectively. 
We propose to show that the solutions for p{”, gf, and r are the same as 
(40) if »=k. 

Consider the set of differential equations obtained by equating the 
coefficients of y* in (11’). The left members are the same as (12A), except 
for the superscript k on p,, q,, and 7,. Let the corresponding right members 
be denoted by P“, Q®, and R®. Then it is found that 


k+1 k+1 k+1 
P® om eta kt} x Q (#) —¢ e ary — eat 


where U{}), V{P, and WS? are functions of the same form as uy, vy, and wy, 
respectively. The complementary functions of these equations are the same 
as at each preceding step of the integration, and the complete solutions may 
also be obtained in the same way as equations (38) were determined. The 
equations analogous to (36) and (37) differ from (36) and (37) respectively 
only in the superscript & instead of 1. The complete solutions will therefore 
have the same form as (38), except that the respective particular integrals 
are functions similar to P®, QQ, and R® instead of P®, Q, and R® as in 
(38). From condition (C,) the arbitrary constants will all be zero, except 
K{, and from condition (C,) this constant is found to carry the factor e*a{t’. 
Hence the desired solutions for p, qg\, and r are the same as (40) if »=k. 


This completes the induction. 


? 

| 
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When these solutions for p{”, gf, and rf, (v=1, 2,....0) are substi- 
tuted in (34) and (33), we obtain as the solutions of (11), 


-k 
| 
j 
(ct) aly’, 


where the superscript on the functions u(t), v@(7t), and w(t) has been 
made to conform with the corresponding power of y. Since the two arbitrary 
parameters a, and y occur only in products as indicated, we may suppress 
either without loss of generality. Let us suppose that a,=—1. 

Equations (41) are the asymptotic solutions of (11) which approach zero 
as t and therefore ¢ approach +o. In the same way there could be obtained 
the asymptotic solutions which approach zero as t and therefore ¢ approach 
—o. It is not necessary to construct these solutions, however, since they can 
be obtained directly from (41), as we proceed to show. 

It was shown in the two paragraphs following equation (20) that if a set 
of solutions of equations (11) is known, then another set can be obtained by 
changing the signs of t, g,, and r, in the former set. On making these changes 
in (41) we thus obtain the asymptotic solutions of (11) which approach zero 
as tT approaches —o, viz., 


2 


7=1 k=1 
3 


j=1 k=1 


Upon substituting (41) and (42) in (10), and returning to the original 
variables £, 7, through the substitutions (4) and (2), the asymptotic solutions 
of the original differential equations (1) in terms of t are found to be 


j= 


4=0 2 ( t) ely’, 


k= 


(43) 


Where the double signs occur, the upper signs give the solutions which 
approach the periodic orbits as t approaches +o, and the lower signs give 
the solutions which approach the same periodic orbits as t approaches —o. 


. 
j=1 k=1 
} 
| 
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Asymptotic Orbits of Class B. 


We consider now the orbits which are asymptotic to the periodic orbits of 
Class B. The method of constructing these orbits is entirely similar to the 
one used in constructing the asymptotic orbits of Class A. 

On putting ;=2 and considering the differential equations obtained by 
equating the coefficients of the various powers of y in (11’), that is, in (11) 
when (33) and (34) have been substituted and y divided out, we obtain solu- 
tions which will give the asymptotic orbits, provided the arbitrary constants 
of integration are so chosen that conditions (C,) and (C,) are satisfied. 

The differential equations arising from the terms independent of y in 
(11’) are the same as (12B), except for the superscript 0 on p,, q,, and f.. 
Their solutions are therefore the same as (31) and (32). When the conditions 
(C,) and (C,) are imposed, we obtain 


pS = ae Pa gS = —a,e 


From the coefficients of y to the first degree in (11’) we obtain the differ- 
ential equations which define p{”, g%, and rf”. These equations have the same 
left members as (12B), except the superscript 1. Let the respeetive right 
members be denoted by and RY’. Then 

= ease "UD, eaze RY? =0, 
where U? and V{% are functions of the same form as wu. and Vz. respectively. 
The complementary functions of the differential equations are the same as 
(31) and (32), but let the arbitrary constants be denoted by J, ...., 19. By 
using the method of the variation of parameters to obtain the complete solu- 
tions we have 


Un) + ( —P2tloo + + IM 

+12 [tig +L + Tl) ] =P”, 
+ IP + L723) =0, | (44) 
err + Vn) (P2U22— Vee ) +19 Veg 
+1 Wo = 0, 


The last two equations are independent of the first four. Since the deter- 
minant of the coefficients of J and J in these last two equations is different 
from zero, viz., A;, the only solutions for J and IM are JO =1=0. 


| 
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The determinant of the coefficients of J, ...., 1 is different from zero, 
viz., A, equation (30), and therefore the first four equations of (44) can be 
solved for if, ...., 1. The resulting solutions are 


(1) 
(45) 


where A® is the determinant obtained by replacing the elements of the k-th 
row of A, with 0, P{, 0, and Q@ respectively. Since the right members P{? 
and Q® do not contain terms in e+’, the integrations of (45) for Jf and 1 
will not contain terms in + explicitly. Such terms arise, however, in the 
integration for / and 1{, but when they are substituted in the complementary 
functions they cancel off. The complete solutions of the differential equations 
at this step are thus found to be 


pP = ey, + LY + + (eg + ugg) + (7), 
= — + LE + (Vag + + eage (7), ¢ (46) 
where L® ...., L@ are the constants of integration, and u(t) and v(t) 
are functions similar to u,, and v,, respectively. From condition (C,) we have 


LP =LP =LP=LP=LP =0, 


and from condition (C,) 
LO (0). 


Then the solutions (46) become 
gs” + (7) ], 
=0, 


where u{)(t) and v{?(t), (kK=1, 2), are similar to and respectively. 
The remaining steps of the integration can be carried on in the same way. 

Proceeding by induction to the general term, we find that 

y+l 

k=1 k=1 
where the functions u(t) and v§?(t) have the same form as wu and vy 
respectively. When these terms are substituted in (34) and (33), we find the 
solutions of (11) to be 


X De aly’, g.= (7) r=0. (47) 
k=1 


j=1 k=1 


As in the solutions (41) we may put a,=1. These equations (47) are the 


| 
| 
| | 
| 
— 
| 
| 
| 
| 
| 
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asymptotic solutions which approach zero as t approaches +. By changing 
the signs of + and q, as in equations (42), we obtain 
j=l k=1 j=l k=l 
which are the asymptotic solutions approaching zero as t approaches —o. 
Since r,=0 in (47) and (48), the asymptotic orbits are of two dimensions and 
are complanar with the periodic orbits of Class B. 
In terms of the original variables &, 7, f, these asymptotic solutions of 
equations (1) are 


&=£,1+%,+ (+7) ely’, 
j=1 k=1 
j=1 


These solutions approach the periodic orbits of Class B as t approaches +0 
or —o according respectively as the upper or lower signs are taken. 


§9. GromETRICAL CONSIDERATIONS. 


The asymptotic solutions (43) and (49) contain the three undetermined 
constants ¢,, «, and y. The constant ¢, denotes the initial time and may be put 
equa] to zero without loss of generality. The parameters « and y are the 
respective scale factors of the periodic and asymptotic orbits. The physical 
interpretation of the parameter ¢ has already been referred to in the latter 
part of §4. From the way in which the initial conditions were chosen in (C,), 
it is evident that ey denotes the £-component of the infinitesimal body’s initial 
displacement from the periodic orbit. 

Let us now consider the directions in which the asymptotic orbits approach 
the periodic orbits. We shall discuss only the orbits which approach the 
periodic orbits of Class A as t approaches +o, since the discussion is 
entirely similar for the other asymptotic orbits. 

As t becomes very large, the most important terms of the solutions of 
(11) are those in e~”’. These terms arise only from the complementary 
functions of the differential equations which define the various p{”, gf, and 
r( in (34). Neglecting the explicit values of the constants of integration 
which are associated with e~"” at each stage of the integration, we find that 
the predominating terms of the solutions as t approaches + are 

r= —e wy [KP yt 
where K®, K®, ...., are the constants of integration. 
14 
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The projection on the &y-plane of the direction of the approach is 


limit 2% — limit pi? 
r=to dt dt OUy— 

This limit is independent of y, but, since u,, and v,, contain sines and cosines 

of multiples of V At, it is indeterminate. The projections on the &¢- and 

7¢-planes are likewise found to be indeterminate. Similar results are obtained 

for the orbits which approach the periodic orbits of Class A as t approaches 

—o, and also for the corresponding orbits of Class B. 


$10. In tusrrative ExaMp.es. 


We shall conclude this article with numerical examples and diagrams of 
the periodic and asymptotic orbits which have been discussed in the preceding 
sections. In these examples the ratio of the finite masses is ten to one, or 
1—u=10/11 and uw=1/11, being the ratio used in the particular periodic 
orbits already mentioned in § 4. 

The solutions for the asymptotic orbits have been carried out to the third 
degree in e for orbits of Class A and to the second degree in ¢e for orbits of 
Class B, but only to the first degree in y for both classes of orbits. 

In the numerical results that are to be found in the tables which follow, 
e has the value 0.5 for Class A and 0.01 for Class B while y=0.1 for both 
classes. The values of ex,;, ey;, and ez; (j=1, 2) in the tables are the 
coordinates for the periodic orbits for the various values of t indicated. The 
values for eyp,, eyg;, and eyr; are the amounts which must be added to the 
coordinates for the periodic orbits in order to obtain the asymptotic orbits. 
For j=2, 

The periodic orbits are represented by the heavy lines in the diagrams, 
and the asymptotic orbits by the dotted lines. The arrows indicate the 
direction of motion. No arrows appear in the periodic orbits in Figs. 2, 5, 
and 8, as in these projections the infinitesimal body oscillates up and down 
along the same curve. The origin of coordinates is taken at the equilibrium 
point marked in the diagram, and the axes are parallel to the rotating &¢ axes 
(see equations (2)). The unit of measurement is indicated in each diagram. 

The following results have been obtained. 


Orbits of Class A. 


Equilibrium point (a). 
ep, cos 2VAtT—1.08 sin 2VAt)+....], 
eq, [—0.75 e+e? (0.22—0.61 cos 2VAt—0.56 sin 2VAt)+....], 
er, [e?(1.58 cos VAtT+0.91 sin VAt)+....], 
pi =1.83—0.23 A=2.548. 
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TABLE 3. 
e=0.5, y=0.1. 

0 — .041 0 0 .0500 .0420 .0400 

in — .043 — .008 .043 .0385 .0369 .0359 

2 — .047 — .015 .085 .0295 .0309 .0315 

— .054 — .021 .125 .0223 .0259 .0270 

4 — .064 — .026 .163 .0173 .0212 .0227 

5 — .074 — .028 .198 .0136 .0173 .0187 

— .084 — .028 .228 .0111 .0142 .0155 

— .095 — .025 254 .0093 .0107 .0118 

8 —.104 —.021 .274 .0080 .0083 .0091 

9 —.111 — .015 . 289 .0071 .0064 .0067 
1.0 —.116 — .007 . 297 .0065 .0048 .0048 
—.115 .009 . 296 .0056 .0029 .0020 
1.4 — .103 .022 | .0045 .0020 .0003 
1.6 — .083 .028 . 223 .0034 .0015 |—.00056 
1.8 — .062 .028 .157 .0025 .0013 |—.00091 
2.0 — .046 .014 .078 .0017 .0011 |—.00095 
2.4 — .048 — .016 — .092 .0005 9 .0006 |—.00065 
2.8 — .087 — .027 — .233 .0002 3 .0002 7 | —.0003 0 
3.2 —.116 — .006 — .298 .0001 4 .0000 9 |} —.0000 9 
3.6 —.102 .022 — .269 .0000 92 .0000 7 | —.0000 4 
4.0 — .060 .024 —.151 .0000 48 .0000 3 | —.0000 18 
4.4 — .041 — .003 .014 .0000 15 .0000 2 | —.0000 16 


The projections of the above orbits on the coordinate planes are given in 


Figs. 1, 2, and 3. 
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Equilibrium point (b). 

ep,=e-""y [e+e (0.222—0.222 cos 2VAt—0.454 sin 2VAt)+....], 

eg, =e"y [0.40 e+e7(0.089—0.182 cos 2VAt—0.037 sin 2VAt)+....], 
er, =e [e(—0.345 sin VAt—0.523 cos VAt)+....], 


pi =3.366—2.837 &+...., A=6.510. 
TABLE 4. 
e=0.5, y=0.1. 

0; .0135 0 0 .0500 .0188 —.0131 

1] .0147 .0022 .0474 | .0300 .0119 — .0094 0 

.2] .0178 .0038 .0917 | .0186 .0073 8 — .0062 6 

3] .0221 .0045 .1355 | .0118 .0051 7 — .0038 6 

4] .0266 .0040 .1602 | .0077 9 .0034 5 — .00226 

.5| .0300 .0025 .1799 | .0051 9 .0022 7 — .0009 57 

6) .0315 .0004 .1878 | .0034 9 .0014 7 — .0005 76 

.7| .03807 |—.0019 .1837 | .0023 1 .0009 36 — .0002 24 

.8] .0278 |—.0036 .1677 | .00149 .0005 74 — .0000 454 

9] .0236 |—.0045 .1594 | .0009 34 .0003 47 + .0000 350 
1.0} .0206 |—.0044 .1047 | .0005 83 .0002 13 + .0000 616 
1.2} .01386 |—.0007 .0150 | .0002 04 .0000 762 |+.0000 497 
1.4} .0166 .0034 | —.0784 | .0000 750 .0000 309 |+.0000 363 
1.6] .0252 .0043 | —.1519 | .0000 308 .0000 136 |+.0000 0936 
1.8} .0312 .0011 | —.1867 | .0000 137 .0000 0590 |+.0000 0257 
2.0} .0289 |—.0032 | —.1837 | .00000594 | .0000 0232 |+ .0000 0056 6 
2.2} .0205 .0044 | —.1169 | .0000 0233 | .0000 0081 0/—.0000 0021 6 
2.4} .0140 .0014 | —.0301 | .0000 0086 6] .0000 0030 8}— .0000 0002 2 
2.6] .0156 .0029 .0645 | .0000 0029 9} .0000 0012 4;— .0000 0001 9 


The projections of the above orbits on the coordinate planes are given in 


Figs. 4, 5, and 6. 
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Equilibrium point (c). 
ep, [e +e? (—0.104-+0.104 cos 2VAt—0.817 sin2 VAt)+....], 
eq, =e"y [3.08 e+e? (—0.320—0.787 cos 2VAt—0.106 sin 2V/At) +....], 
er, [e(—0.719 sin VAt—3.145 cos VAt)+....], 
=0.457—0.0252+...., A=1.082. 
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TABLE 5. 
e=0.5, y=0.1. 

0 .0050 0 0 .0500 . 1402 — .0786 
2 .0099 .0228 .0989 .0418 .1279 — .0733 
4 .0238 .0418 . 1939 .0349 .1180 | —.0657 
.0443 .0536 2803 .0297 . 1098 — .0560 
.8 .0678 .0563 .0264 .1030 — .0455 

1.0 .0905 .0493 .4138 .0245 .0960 | —.0346 
1.2 .1084 .0340 .4550 .0236 .0893 | —.0239 
1.4 .1185 .0129 .4766 .0234 .0825 — .0139 
1.6 .1190 .0105 .4781 .0233 .0754 | —.0051 
1.8 .1099 .0320 .4584 .0229 .0680 .0026 
2.0 .0927 .0481 .4190 .0221 .0612 .0088 
2.2 .0703 .0560 .3619 .0209 .0543 .0135 
2.4 .0466 .0543 . 2890 .0191 .0481 .0168 
2.6 .0256 .0434 . 2035 .0169 .0426 .0187 
2.8 .0110 .0251 .1094 .0147 .0377 .0194 
3.0 .0051 .0025 .0105 .0123 .0341 .0191 
3.4 .0220 .0400 — .1843 .0086 .0289 .0161 
3.8 .0653 .0564 | —.3475 .0064 5 .0251 .0113 
4.2 . 1069 .0359 — .4517 .0057 3 -,0218 . 6060 
4.6 .1194 .0080 — .4790 .0056 2 .0184 .0014 
9) .1101 .0318 — .47238 .0051 0 .0151 .0015 
0.4 .0491 .0549 —.2971 .0046 6 .0119 — .0040 
5.8 .0121 .0272 | —.1195 .0035 9 .0093 | —.0047 
6.2 .0081 .0182 .0787 .0014 5 .0076 | —.0044 


The projections of the above orbits on the coordinate planes are given in 


Figs. 7, 8, and 9. 
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Orbits of Class B. 


Equilibrium point (a). 
[e+e (—1.20+1.20 cos ot+3.95 sinot)+....], 
[0.75 (0.90 +4.65 cos ot +2.46 sin ot) +....], 


=1.833+()e+...., c=1.68. 
TABLE 6. 
e=0.01, y=0.1. 
0 .0100 0 .0010 00 .0007 13 

2 .0094 — .0088 .0007 00 .0004 77 

4 .0078 — .0166 .0004 90 .0003 31 

6 .0053 — .0225 .0003 42 .0002 31 

8 .0023 — .0259 .0002 38 .0001 63 
1.0 .0011 — .0264 .0001 63 .0001 16 
| .0043 — .0240 .0001 13 .0000 820 
1.4 .0070 —.0189 .0000 776 .0000 581 
1.6 .0090 —.0117 .0000 534 .0000 412 
1.8 .0099 — .0031 .0000 361 .0000 289 
2.0 .0098 .0058 .0000 248 .0000 203 
2.2 .0085 .0140 .0000 169 .0000 140 
2.4 .0063 .0207 .0000 117 .0000 0970 
2.6 .0034 .0250 .0000 0809 .0000 0664 
2.8 .0001 .0266 .0000 0560 .0000 0452 
3.0 .0032 .0252 .0000 0392 .0000 0306 
3.2 .0062 .0210 .0000 0275 .0000 0208 
3.4 .0084 .0144 .0000 0193 .0000 0141 
3.6 .0097 .0062 .0000 0136 .0000 0095 
3.8 .0100 .0027 .0000 0094 .0000 0065 


These orbits are shown in Fig. 10. 
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Equilibrium point (b). 
ep, =e [e+ e(0.69—0.69 cos ot—2.47 sinot)+....], 
=e [0.40 e+e?(0.28+4.60 cos ot+3.07 sinot)+.... 


Equilibrium Points in the Problem of Three Bodies. 


o=2.61. 
TABLE 7. 
e=0.01, y=0.1. 
0 .0100 0 .0010 00 .0004 49 
wa .0097 .0103 .0007 10 .0003 70 
a .0087 .0199 .0005 04 .0002 34 
3 .0071 .0282 .0003 59 .0001 66 
4 .0050 .0345 .0002 55 .0001 17 
.o .0026 .0385 .0001 82 .0000 83 
€ .0000 4 .0399 .0001 31 .0000 58 
7 — .0025 .0386 .0000 933 .0000 398 
8 — .0050 .0347 .0000 670 .0000 275 
9 — .0070 .0284 .0000 481 .0000 190 
1.0 — .0086 .0202 .0000 345 .0000 131 
1.2 — .0100 .0003 .0000 178 .0000 063 
1.4 — .0087 .0196 .0000 092 .0000 031 
1.6 — .0051 .0344 .0000 047 .0000 016 
1.8 — .0001 .0399 .0000 024 .0000 0087 
2.0 + .0049 .0348 .0000 012 .0000 0047 
2.2 +- .0086 .0204 .0000 006 .0000 0026 
2.4 .0100 .0006 .0000 003 .0000 0014 


These orbits are shown in Fig. 11. 
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Equilibrium point (c). 
(25—25 cos cot —23 sinot)+....], 
[3.08 e+ —56 cos ot +31 sinot)+....], 


p2=0.475+()P+...., o=1.07. 
TABLE 8. 
e=0.01, y=0.1. 

0 .0100 0 .0010 0 .0032 9 

e .0099 — .0022 .0009 3 .0031 7 

2 .0098 — .0043 .0008 7 .0030 6 

3 .0095 — .0063 .0008 1 .0029 7 

4 .0091 — .0083 .0007 6 .0028 7 

.o .0086 — .0103 .0007 2 .0027 9 

6 .0080 — .0120 .0006 8 .0027 0 

i .0073 — .0137 .0006 5 .0026 2 

.8 . 0066 — .0152 .0006 3 .0025 4 

9 .0057 — .0165 .0006 0 .0024 6 
1.0 .0048 — .0176 .0005 8 .0023 9 
1.5 — .0003 — .0199 .0005 4 .0020 5 
2.0 — .0054 — .0169 .0004 6 .0017 1 
2.5] —.0089 — .0090 .0004 2 .0013 7 
3.0 — .0100 .0014 .0003 7 .0010 6 
3.5 | —.0082 .0140 .0003 0 .0007 8 
4.0 — .0042 .0183 .0002 3 .0005 7 
4.5 .0010 .0200 .0001 4 .0004 1 
5.0 .0060 .0162 .0001 2 .0003 0 
5.5 .0092 .0078 .0000 8 .0002 3 
6.0 .0099 — .0027 .0000 6 .0001 9 


These orbits are shown in Fig. 12. 
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Concerning the Invariant Theory of Involutions of Conics. 


By Wayne SENSENIG. 


§1. Lytropvuction. 

The complete system of invariants of a single conic 
consists of four members: f=a?, L=a2, D=a?, and u,. The complete system 
of invariants of two conics consists of twenty forms and this was derived first 
by Gordan and was first published in Clebsch’s “ Vorlesungen tiber Geometrie.”’* 
H. F. Baker derived the complete system for three, and incidentally two, 
conics.t Ciamberlini has published a determination of a complete system for 
three conics. { 

The object of this paper is to derive and reduce in terms of the system 
of two conics a2, b2 the complete simultaneous system of the involution 
K=a?+kb? and the harmonic conic F= (a6x)*. The complete system of K and 
F is expressed in terms of the complete system of the two base conics f=a? and 
g=b?, which determine the configuration. 

In Bochers “ Introduction to Higher Algebra,” p. 164, it is shown that if 
two conics f, g intersect in four distinct points there exists a (non-singular) 
collineation which will reduce f and g simultaneously into the normal forms 
given below. Note that the equalities hold only by virtue of the equations of 
the collineation : 

(1) 
In Table I, I give the well-known symbolical forms of the system of f and g, and 
also the normal system which I have computed from these symbolical forms by 
particularizing the coefficients of f and g, respectively, according to the follow- 


ing scheme: 


Throughout the paper we have used the algorism for ternary transvection 
given by O. E. Glenn in the Transactions of the American Mathematical Society, 


*Cf. Osgood, AMERICAN JOURNAL OF MATHEMATICS, Vol. XIV (1892), p. 262, “System of two 
Simultaneous Ternary Quadratic Forms.” 

+ Cambridge Philosophical Transactions, Vol. XV (1889-93), p. 62. 

t Giornale di Matematiche, Vol. XXIV (1885-86), p. 141. 
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Vol. XVII (1916).* This algorism is analytical but it will suffice here to 


describe the transvectant t of index te of the two forms 


where A and B are constants, as AB/u times the sum of all terms gotten by 
forming in the product $1, in all possible (u) ways, i convolutions of the type 
Qiz,,j of type (a,b,u), k of the type ay,,, and J of type eae Then 7 is 
abbreviated as 
t= (9, 

In deriving the system of K and F' we have evaluated each member as a 
transvectant and then reduced each form to a rational integral function of the 
members of the fundamental system of a? and b?, expressing each one in the 
form of a polynomial in k. These reductions were performed by means of 
the following fundamental ternary identities, and their modifications: 


(abc) d,—(abd)c,+ (acd) b,— (bed) a,=0, | 
a, 60, wu 
a a a 
(abu) = ag bs Ug ( 


The result for each case was then checked by direct verification, for the system 
of normal forms given in Table I. Table II contains the summary of the results 
of the reductions for all twenty concomitants in the system for K and F. 

In Table III of the paper we give the fundamental system of the involution 
a?+kb? taken with a third general conic c?. All concomitants in this system 
are expressed as rational integral functions of certain forms of the system for 
three conics as given by Baker; in fact in terms of the sixty-one forms of this 
system which are the simplest and therefore the most important in geometrical 
applications. This set of sixty-one concomitants is given below: 


f=az, g=b:, ®,=(acu)’, P=(beu)’, Cy=(acu)a,c,, Cy = (beu) b,c, , 
‘OF =a,b,(abu), Ag=t, A,=¢, Cy 5 
F,=(ayx)’, F,=(By«)’, (aye), 
Cy=0,0,0, (aya), L'=y5, 
T,=b,cgugu,(beu), J=u,, M=(abc) (abu)c,, 

N =(abc)*, 0,=a,b,y,(abu), O,=a,b,y,(abu), O,;=b,c,a,(bcu), 

0, (acu), P=a,a,b,b,, Q,=4,b,c,y, (abc), (abc), 


* Cf. Glenn, “Theory of Invariants ” (Boston, Ginn & Co., 1915), p. 227. 
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R, =(a’be) (acu) (abu)a,, R,=(ab’c) (bcu) (abu)b,, S,=b,c,(beu), 

8, =agCg(acu), S,;=a,b,(abu), T,=b,b,a,(ayx), T,=bbay,(ayx), 

T; =a,0,8,(Byx), T,=aga,u,(Byx), Ts=b,b,8,(Byx), 
U,=aga,(Byx), V=a,b,c,(abc), W=a,b,¢,(abe), 

X= GC p4,C,(Byx), Y=(acu)(abu)(beu), Z,=a,cgugu,(acu) ; 
where those preceding J are concomitants of two conics and those following 
J, of three conics. 


TABLE I. 
Normal System of Two Conics. 
g 
D =a=Gabe. 


L =a2=2 (bcu?+acuz+abuz). 
L' 
= (az, (abu)?= (br+cq)ui+ (ar+cp)ug+ (aq+bp) us. 
C, =(az, bz) (abu) a,b, 
= (aq—bp) (Cp—ar) (br—cq) . 
9. Aye= (az, (agr+bpr+cpq). 
10. C, Bi) (agra, + 
11. Cy, =(a2, (bepxu,+ 
12. Aye= (ar, 
13. Cy, =(a1, (aBx) 
=4 [cr (bp—aq) bg (ar—cp) + ap (cq—br) ugis%]. 
14. F Bi)o= (aba)? 
=4[ap(br+cq)aj+bq(ar+cp) cr (bp+aq) x5]. 
15. Cs = (a2, 
=4[agr (br—cq) xu, + bpr (cp—ar) cpg (aqg—bp) 
16. Cy = (a2, 
=2[p’(cq—br) + 9° (ar—cp) (bp—ag) 
17. C, (abu) a,a, 
18. Cs =(a%a2, 2) 
=4[abr(aq—bp) acg (cp—ar) + bep (br—cq) 
19. =(fL, 
=4[a*qr (br—cq) +b*pr (cp—ar) +c*pg(aq—bp) ] 
20. T =(fL, 048, 
=4[p*be(cq—br) +q’ac(ar—cp) +1’ab (bp—agq) 
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§ 2. Repuction oF THE System or K anp F. 


The complete system of K and F is obtained by forming and reducing the 
transvectants of K and F which correspond respectively to the transvectants 
of f, g required to produce the system of f,g (cf. Table I). Thus a member 
of the system of f, g is ®=(f,g)%. Hence the corresponding form of the 


required system is 
@ = (K, F)§. 


I have not included in this paper all of the reductions which were necessary 
to produce Table II; but only typical ones. The most complicated cases of all 
were C;, G’, and I’. Weadd that the reductions for 3,4,5,6 in Table II had 
been given previously.* 

The Form ©’. 
(K, F) (a (a8)2)0= (42, (a8) 2) 00 
= (aaBu)?+k(baBu)?= (agug—agua)* +k 
= DL’ + (aa’a’”) (a’a’'u) agus 
+k [D'L + (bb'b") (b’b’'u) 


The Form C;. 


(a2 + 2k (yd)2)%, where a=y and 
= (ah, (v8) 2) 2k[ (ab), (yd) (v5) 
= (ay8) a4 (8x) +2k (aby) (abu) 
1. (ay3)a4(y8x) (wyS) {(y82)ay— (adx) =4 [(ayd) 
(ay8)*2,—# (aya). 
aj 
= — )] ab’ — (a’a" a" Jay’ 16, 
+ ag’ at" 15, 
=D - - al’al’B,+D - - 


*«“ Algebra of Invariants,” by Grace and Young, p. 293, Cambridge University Press, 1903. 
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2. (abyd) (abu) (8x) = (abu) (a Bx) =a,bs (abu) (a8) 
—agb,(abu) (a 8x) =a,bg(abu) —k, 
=4DD'u,—hy, since a,b,(abu) =agb,(abu) 
= (abu) (abu) = (aa’a’”’) ayy) bgb; 
— (aa’a’’) (a’a’'b’) bgby =4D bu) (bb’b’’) b; 
—4D (b’bu) (bb’b”) by =4D (bb’'u) (bbb’) + 4D (bb’u) (bb’b’’) by” 
=$DD'u,+3DD'u,=$DD'u, 
38. (By8) (yea) B.=4 (By8) [ B,+ (Byx)8,] (By8) [ (Byx) 
=4 (By9) [ (By8) + (88x) y.] =4 (By8)*a, +4 (Bys) (Bdz)y,. 
=$ [a gD, pb a, + 01, | 
= a’ (yng y= a") bs a, 
= a, = 
(ByS) By kg=agb, (abu) (aBx) DD'u, 
+ 112A 190, + 
= DC, + 2k DD'u,—Ay2C 2+ Anz 3} 
+h? 


The Form C}. 


1. fata®, (a’B)2|% = (aa’B) (aa’Bu)a,a, = (aa’B) 
=— (aa’B) 3,a,0,=— (aa’’’a'’) Ba) a,a,6, 
=—43D (Bax) a,8,=3DC,. 

2. (aB)2}% = (aaBu)a, (abu) 
= (aga, a, bu) 

(Terms containing a,’b, vanish.) 
For, —ag” a,b,a, (abu) 
=a, 'a,b,0,} (a’’bu) ag— (aa’’b) ug— (aa’’u) bg] 
=4b,a,ug(aa’’b) — (bb’b”) (aa’”’u) 
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But, * (xB) (4,0) a,— a5] 


= b,] aya, 


)agd, tig —4D (ba’u)b, Ag Ug=$D (abu) agb,ug=3DC,. 
3. {a262, (aaB'u)a,8,= (Baf") | 4,8, 


4. {bat, = (aa’B) (ba’Bu)b,a,= (aa’B) 


=b,ga,b,a, = %DC,. 


5. {b2(ab’”)2, (a3)2|% = (aba) (babu) b, (ab’”u) 
= (bga,—b,68,) b, (ab’’u) (ab’”’u) 
+ (ab’’’u) * = —34D’C,+ 
=Ay,.C,—4A,.C,; * similar to 2. 
6. (aB’)2}% = (baB’u)b,B,= (bby 
= (BBY) Br (bub’)b/b,a, 
=4D' (bub’) a, [b{b,—b/b,] =4D' (b’bu) a, (b’bax) 
Also, by b’~b", and 
Hence 
+k? (—4D'C,). 
The Form G’. 


(a, +2k (ab), +h Bt), (a (ab) 

+h (dian, i+ 2k (bz 2, 
+ (a8')2 (ab) 
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(aa’B) (a8'x) a, (a’Bx) + a, (a Bar) 
+ ay, (882) a,(aB'x) +kby (aa’B) (aB'x)b,(a' Ba) 
+214» (ab™aB) (ab”’B'x)b, (aBx) 
+3b3b2 [dy (aa’B) (aaa) a, (a’Bx) +2ka, (a’”baB) a, 
+ Kay (Ba'x) a, + (aa’B) (aaa) b, (a’ Bx) 
+2k?by (ab’aB) b, (a Bax) (BaB') (Ba’x)b, 
—2a2b3[ (aa"b) (aa’B) a, (a Ba) 
42k (aba) a, (Bx) 
+l? (aa’"b) (Ba’’ba) a, +k (bab’”)(aa’B) a’ Ba) 
42k? (ab’"aB) b, (a Bar) 
+-k*(bab”) (Ba(3") (Bab¥x)b, (a8"x) 
=3D +h?(—3D'G) +k (0) +24 (0) +°(0)] 
(0) +2(0) +h(3D@) +1 (0)] 
—2DD’[—4DG+ (0)+% (0) 


(8) 
Collecting terms, 


G' = §3D°D' +k (4DD' Ay.—6D +h? (6D' Aj,—4DD' Ay.) —3 DD" - |G. 
The calculations are shown below: 


1. ay (aa’B) a, (a’Bx) [ — 
=4ay0,(aa’B) (aa’x) (8’Bx) =4a,(aa'x) [ (aa’B) Bi,—(aa’B’) By] 
=4a,(aa’x) (88'x) [ (a’BB’)a.— =0. 

2. by (aa’B) b, (a’ Ba) = za) b, (aa’B)(a’ Bx) 

=4D’ (xaxz), ete. =0. 


3. Aq (aa’B) (aa’’x) a,(a’ Bx) = etc. =0. 


4. Dan (aa’B) (00%) b, (a’ Bx) =Dan(aa’B) 
=ban(aa’B) b, Bx) (aa’B) azay.b,(a' Bx) 
= [ayag Bx) 
— (a%a"'b) [ayag —agay,] a,b, (a’ Bx) Bx) 
Bx) —4Db,a,a,b, Bx) Bx) 


Bx) (bb’b”) (b'b" xa’) %DG—0. 


i 
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5. 


10. 


(aa°b) (aa’B) a, (a’Bx) 

— (a' Bx) a,(aa"b) b,(a’ Bx) a,(aa‘b) 

( Ba) 4g 0,0,(aa"b) —0 

=—4D(azay +4D agb,a,(aab) 

=4Dazazb,(aa‘b) {a,a3—agay | aza,—ayag} 

—4Dazagb,b, =—4DG—4{DG=—4DG4. 


(bab’”) (aa’B) b, 
= (a’Ba) (bab””) 
= (a’ Bx) b, (bab’’) (a’ Bx) b,(bab’’’) +0 
aa,b,(bab”) = 3D 
= — = —3D (b'"b’b”) (b'b"a") a,b, =0. 


— 
dg (a’’baB) (a’"bB'x) a,(aBx) bq) (ag b,—az by) (a 8x) 
= (a’’’a’a’’) ( a’a’” Bx) Agr (a’a" Bx) 
— b,Aga,(a8x) + (b'"ba) a, (482) 


by (ab’”B’a) b, (aBx) = , ete. = 0. 


dy (a’"ba'x) a, (48x) = (aa’"a"”) a, 
ete. =0. 


But, by aga, baby’ ] 
[ (Bar) — ( a’ ayn] 
=42,b,d, [ (aBa’) — (axa’) =0. 
b, (a Br) = 
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11. (a’"baB) a, (a Bx) 
For, (aBx) a,j b,— az’ } 
(aa'’b’”’) pa, )b,=0 
and (aa'’b’”’) a,’b,a,(aBx) (ba'’b’”’) a,’ =0. 
But, — (aa'’b’”’) b, 
= [ (aal’b’”) ag’ — 
=—4b,a,(aBx) [ bg’ + ag] 
| a,b, bbs’ (a8x) 
=—44y.G. Also a,’ 
(aBx) { (aat’b’’) b,— (aal’b) by} 
(a Bx) { (abb’”’) — ag 
(bab) (abb’”) +4 (a2) 
a’ Ba) a, a, '(abb’’”’)= 4D a,’ (abb’’”’) 
(abb’”’) 
DD’ DD’ (a’"b'"a) =0. 
12. (ab"aB) (ab b, (aBz) 
+ agby” (aBx) a,b, 
= (aBx) (aa’’’b"’) a,— (bab!) ay’ } 
+$a,(aBx) a,b,{ by” — (ba’"b’”) BEY} 
+44,(aB2) (b'"a'"b™) a,b, | (bbb) ag’ bY) bg} 
=4 AinG. 
13. ay(Bap") a, (aB'x) | yn} (B'aa)a, 
by, 
14. (BaB") b, (aB'x) = (Bap) =0. 
15. ay(BaB") = (aa’"a'”) =0. 
16 


. 

1! 

| | 
| 
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16. by (BaB’) (Ba’x)b, (aB'%) — (B’xa) b, 
—4D' (xaa’) = D’ (aa’x) 


(bb’a’)b,b,a,=0. 


17. (aa’”’b) (Baf") a, 
= (aa’'"b) (a8'x) | 
= (aa’’b) (a8'x) a, + | 
— ba, (aa’’d) 
— (b’b'"b'") bya, (aa’’b) 
=4D’ (bb’xa) (aa’"b) —4D’ (a’"b'xa) 
—4D' (aa’"b) +4 bia, (aa’’b) 
= —¥ 
=% D’ (a’"a'a"”) (a'a’"b’) 
= (b’ba) bib,a,=0. 


18. (bab) (Gap) 
= (bab") (Bab¥ax) b, 
= (bab’) b, —D Aga, +b (a,b3—agby ) 
= a’) Doren a, (bab) b, 

— a’) by (bab’)b,, ete. 
=4D'(a’'b‘b’) ) —4D’ (a’ab’) (bab’), ete. 
=} [ (bab’) bj — (bab’) by] —0, ete. 
=} [ (ab’b’)b,— (bb’b’) a,], ete. 
(abYb’)b,, ete. ete. 
], etc. ete. 


Combining these results as in (3) we obtain G’. 


Le 


| 

| 

| 

| 


go 


15. 


16. 


17. 


18. 


19. 


20. 
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K 
F 


. (D)’ 


o’ 
C; 
C; 
Cs 
Ans 


F’ 


G’ 


ad 


(D’)’ 


TABLE II. 
System of a?+kb?, (aBx)?. 
=ai+kbd?. 
= (aBx)’. 
=ai+2k(abu)?+k 


(abu)*}. 
=D+3kAy.+3h Ayo 
DD’. 
=4$DL'+ 
=C,—kC;. 
=$ DD’ Aye) 
=, 
= + 2k ($DD'u,—Ay2C,+ — A 
+h? D'C;). 
=§{DAy,+2k(4 DD! + +h?-$D' Ay. 
}. 
=${3DAj_F +2DD'Ayf—¥ D*D’g 
+2k(3DAing—DD'A of —Ay2DD'9 +3Aj.D'f +DD'F) 
+2h Ay,C;—3D A 
+I? DD"C,) 
=$ Ay»C;—DD' D°D'C, 
+k + DD’ DD"C;) }. 
+h? (—Aj2C,+ D’C,) (—4D'C,). 
=$13D GD + 2DD’Cy—3D A 
+k? 
+k (2 DD?C,—3D' 
D°D’ +k (4DD’Ay,—6D 
=$ {3 D°D'+k (4DD’ Ay,—6D A 


| 
| 
| | 
= | 

mC; 

m C; 

= 

| 
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§3. System or a3+hb?, 
TABLE III. 
(a2-+-kb?2, c2)ii in terms of the system of Baker. 


= , cz) =P, +k. 
= (a2 =CrthCy. 
= (ai +kb:, 00 = Cu 
= [a2+2k(ab)i +h Bi, 
[an +2k(ab)i+h Bi, =As+2kN+hA,. 
Cy =[ai+2k(ab)i +h Bi, 
F” =[ai+2k(ab)i+ Bi, 
+203) +h? (2R,—S.J +20,) 
= Cog +h 
+h? (—24,C,+298,+3T,—T,—JU,) 


=[ (ai + 2k(ab)i+h Bi), 
=G,+k[24,V—2f/W+X,]) [— +k*G,. 
= [ (a3+kb?) (ai+2k(ab)i 10 


+8,L"—28,8,4+2,] +#[— —24,¥ —S,L" + 28;0,+23] +T, e 


| 
| 

| 

| 

| 

| 
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Note on Seminvariants of Systems of Partial Differential 
Equations. 


By A. L. NEtson. 


1. Introduction. 


In the discussion of the projective differential geometry of a geometric 
configuration by means of Wilczynski’s method, one of the necessary steps is 
the construction of a fundamental set of seminvariants, that is, such a set that 
any seminvariant whatever may be expressed in terms of them and their 
derivatives. In the cases where the completely integrable system of partial 
differential equations employed has one dependent variable, this construction 
has been accomplished by the reduction of the system of equations to its 
canonical form, the independent coefficients of which form are the fundamental 
set required. Let us illustrate the process by the case of plane nets.* The 
completely integrable system of equations for this case is the following: 


by, Hey, (1) 


The transformation (2) 
yields a new system of equations of the same form as (1), with the coefficients 


The integrability conditions enable us to find a function p(u, v), such that 
P,=a+b’, p,=a’+b”. We find also that under the transformation (2) 


these combinations become 


* Cf. E. J. Wilczynski, “One-Parameter Families and Nets of Plane Curves,” Trans. Amer. Math. 
Soce., Vol. XII (1911), pp. 473-510. 


Pu=Pu—3 P=P, 
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If, therefore, 4 be so chosen that 


we get the special unique form of (1), which shall be indicated by capital letter 
coefficients, characterized by the relations A+B’=0, A’+B”’=0. This is the 
canonical form of (1). We list its coefficients for later comparison. 


A =}( a—2b’), B C =c +$#a®? +4ab’+4a0'b +4bbd” | 
—ta, 

A’ =3(2a’— 6b”), B’ =4(2b’— a), C’ (4) 
—tab”—3a,—4b,, 
A" =a", B’=}( b’’—2a’), C’=c" + $b’? + + +400” 
—ja,—ta". | 


Only seven of these coefficients are independent. They are the fundamental 
seminvariants. For, consider any system of equations of the form (1), to which 
we shall refer as (s), and any other system (¢), obtained from it by any trans- 
formation of the type (2). Let (s) be reduced to its canonical form (8), and 
(t) to its canonical form (7). The coefficients of (S) are the same combinations 
of the coefficients of (s) as the corresponding coefficients of (7') are of those of 
(t). But, since the canonical form is unique, the corresponding coefficients of 
(S) and (7) are identical. Hence they are seminvariants. Moreover, any 
seminvariant whatever, J, is a function of the original coefficients of (s), 
i. [=I (a, b, ...., a,,....). Since it is a seminvariant, it must be 
identical with the same function of the corresponding coefficients of any system 
of equations obtained from (s) by a transformation (2). In particular, 
I=I (A, B,....,C”, A,, A,, ....). Therefore J is a function of the seven 
independent coefficients of the canonical form and of their derivatives. Hence 
these coefficients are a fundamental set of seminvariants. 

The determination of this canonical form depends upon the ability to find 
such a function py. Green* has shown that this can be done for completely 
integrable systems with one dependent variable and nm independent variables, 
provided certain conditions are fulfilled. However, the expressions for p, and 
p,, which, with their derivatives, must occur in the coefficients of the canonical 


*G. M. Green, “The Linear Dependence of Functions of Several Variables, and Completely 
Integrable Systems of Homogeneous Linear Partial Differential Equations,” Trans. Amer. Math. Soc., 


Vol. XVII (1916), pp. 483-516. 


, 

j 
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form, are not always simple, and as a result the seminvariants which come from 
this regulation method are sometimes rather complicated when expressed in 
terms of the original coefficients. We proceed to indicate how (in most cases, 
at least) an alternative set of simpler seminvariants may be computed. 


2. The Effect of Transformation (2) in the General Case. 


We assume that the completely integrable system of equations (corre- 
sponding to (1)), to which we shall refer as (1’), has one dependent variable, 
y," and n independent variables, u,,....,u,, (%>1). It will then consist of 
a certain number, p, of equations which express certain p derivatives of y in 
terms of g chosen derivatives (which Green has called the primary derivatives). 
Let us also assume that among the q primary derivatives are none of order 
higher than the order of the lowest of the p derivatives in the left members of 
the p equations. 

The transformation (2) gives rise to the following expressions of the 
various derivatives of y in terms of the derivatives of 7. 


(i=1, 2, n), 
D1 2005 
where 


mre = 
= p;! (l,—p,) ! and Our... 


If we substitute (2’) in the system (1’), and arrange the resulting equations 
properly, we get a new system of the same form as (1’), in which the dependent 
variable is 7. The coefficients of this new system, which, in collected form, we 
shall refer to as (3’), are equal to the corresponding coefficients of (1’) plus 
linear combinations of fractions whose numerators are derivatives of A, and 
whose common denominator is A. Each coefficient (3’) receives additional 
terms from those y-derivatives in its equation which come by differentiation 
from the y-derivative belonging to this coefficient. (Thus, for example, in (3); 


b’, the coefficient of 7,, can receive only one additional term, — =, which comes 


from y,,; C’ receives one A-fraction from each of the derivatives y,, y,, Yuy-) 


*This assumption is merely for the sake of simplifying the discussion. The process may be readily 
extended to cover the cases where the number of dependent variables is greater than one. 


i] 

; 
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We arrange these coefficients (3’) in classes as follows: In the zero-th 
class shall be placed those coefficients in each row which are coefficients of 
those primary derivatives having the same order as the derivative in the left 
member of the equation belonging to that row. The 1-th class shall consist of 
the coefficients of those primary derivatives in each row of order 1 less than 
the derivative in the left member of the equation belonging to that row. (Thus, 
in table (3), the coefficients in the first and second columns are of the first class; 
those in the third column are of the second class.) The zero-th class coeffi- 
cients, if present, must be seminvariants, since their y-derivatives can not 
yield by differentiation any other y-derivative in their equations. The first 
class coefficients will be altered, if at all, by the addition of multiples of 
a-fractions of the first order, the coefficients of which A-fractions are integers 
or integral multiples of zero-th class coefficients (3’). In general, the i-th 
class coefficients will be altered by the addition of linear combinations of 
a-fractions of order not higher than 7. Each one of the coefficients of these 
a-fractions will be an integer or an integral multiple of some one coefficient 


(3’) of classes 0 to i—1. 


3. Pseudo-Canonical Forms of (1’). 


Now it very frequently happens that under suitable restrictions we may 
make n first class coefficients (3’) vanish by a proper choice of A. For example, 
in (3), under the assumption a,=b,, we might have made @’=b’=0. Or we 
might have made G=b”=0, provided a,=bi’. If we can do this, under the 
necessary restriction the new form of the completely integrable system deter- 
mined by this choice of 4 is unique, because of the nature of the equations in A. 
We shall call this unique form a pseudo-canonical form of (1'). Moreover, 
the necessary restrictions are always seminvariantive ones. For, consider 
certain such first class coefficients, 


a;,— k; log (é<x1, 2, n), (5) 
where the k; are at worst integral multiples of seminvariant coefficients of (1’). 
In order to choose A to satisfy the equations 


a,— loga=0, (i=1,2,....,m), (6) 


the assumptions necessary are 


(a\ 9 /a 


which are easily seen to be seminvariantive from (5). 


| 
| 
| 
| 
| 
| 
| 
| 
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Let us assume the necessary restrictions (7) and choose 4 in accordance 
with (6). However, there is a precaution to be observed in substituting for A. 
Whenever a cross-derivative, 


Our" (at least two p’s $0), 
1 eee ‘n 


occurs, equations (7) afford us a choice in the substitution. We must, 
throughout a coefficient (3’), replace such cross-derivatives by the proper 
derivatives of the same members of equations (7). The reason for this pre- 
caution will become apparent in the next section. Then, since the restrictions 
are seminvariantive, and since the pseudo-canonical form is unique, the 
coeflicients of this form are seminvariants, by precisely the same argument as 
was made in Section 1 for the coefficients of the canonical form in the case of 
plane nets. 


4. The Coefficients of a Pseudo-Canonical Form are a Fundamental 
Set of Seminvariants. 


Now, if the restrictions (7) be removed, the pseudo-canonical form will 
no longer exist, but the combinations we have obtained are still seminvariants. 
For, let us form the same combinations of the coefficients (3’). If we do this 
under the assumption (7), of course all traces of 4 will disappear. That this 
must also be true without these assumptions is evident from the following 
considerations: 

Let any one of these combinations be denoted by f(a, b,....). It is 
composed of coefficients of (1’) of various classes. However, from the nature 


of the choice of A, the only derivatives which appear in f(a, b,....) are 
derivatives of coefficients of the first class. Now, let us form the same com- 
bination of the coefficients (3’), f(@,b,....), and examine the resulting 


expression. It may acquire derivatives of coefficients of (1’) only from the 
following sources: 1. Those derivatives which appeared in f(a, b,....) will 
remain. These, as we have noted, are derivatives of coefficients of (1’) of the 
first class only. 2. Those arguments, @, b,...., which are of the first class 
(whose expressions are given by (3’)), may have as coefficients of A-fractions 
zero-th class coefficients of (1’). Hence, when the expressions (3’) for these 
arguments are substituted in {(@, b,....), the resulting expression may con- 
tain derivatives of these zero-th class coefficients. It thus appears that the 
only derivatives of coefficients of (1’) of any class except the zero-th which 


17 


| 
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may occur are those from the first source. Moreover, we have taken the pre- 
caution that among these derivatives in f(a, b,....) shall appear only one 
member of any of equations (7) and its derivatives. Hence, if we collect the 
expression resulting from substitution of the values (3’) for the arguments in 
f(a, b, ....), the collection being with respect to the A-fractions of various 
orders, no one of these fractions may have as coefficient an expression which 
vanishes by virtue of assumptions (7). Therefore, since in the resulting 
expression all terms in 4 cancelled under assumptions (7), the same thing 
must happen without these assumptions. 

It is to be observed that in a pseudo-canonical form, there are pg—n non- 
vanishing coefficients, precisely the number of independent coefficients in 
Green’s generalized canonical form. Moreover, these coefficients are inde- 
pendent, since no two arise from the same coefficient of the original system 
(1’). From this it may be argued that the seminvariant coefficients of a 
pseudo-canonical form are a fundamental set. 


5.* Seminvariants Arising from Cross-Derivatives of log i. 


In certain cases there is a lack of definiteness in the formation of semin- 
variants by the process outlined above. As mentioned previously, whenever a 
a-fraction involving a cross-derivative of A occurs in a coefficient (3’), either 
member of the proper equation of (7) may be used in forming the seminvariant 
corresponding to this coefficient. Indeed, certain expressions of the type 


(4; (4; 
m+n 


may be used instead of or 


A coefficient of the type in question may be.put in the form 
Jed, n 
i<i i,j 
+ > a”? D®?+sum of terms of type aD” ,(hy, 


2 
where D@ P=; 2 log ~. Let us denote one of the possible substitutions for 


D“? by [D“?], and the other by {|D“”!. We wish to find what expressions 
of the form 


Fo P= 


(m? (8) 


* Section 5 was added at Professor Wilczynski’s suggestion. 
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may be used to replace D“”. We shall find it convenient to use the notation 


(i, j) (i, 7) (i, 9) (i, 9) (i, 9) i, 


ws (9) 
and it is easily verified that 


Now, if we denote by 6=6+86 the same combination of the transformed 
coefficients (3’) as @ is of the original coefficients of (1’), we see from (5) that 
(D@?] = =| Tf, therefore, F“” is to 
replace D“”, it must happen that 


i) 
(i) 
In order to compute OF ts ivy let us successively differentiate the 
equations 


which result from (8) and (9). We obtain 


Pi=0, 
eee Pn u 1 n 
Pi =0, 


From these equations 


(i, 1 Pn=0, 
l 


n— 
Uy m Dt 
by 


Pe-1 Peri Uy cove 


4, 


Levee Uy Uy 


and will have a similar expression with m””, n®?, D®” and 


replacing m“?, M%?, and respectively. Use of the 
relation 


(i, (i, i) 
OF, 
U1 Un Uy Uy 
| 
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gives us, after considerable labor, 


(i, j) (i, 9) 
oF 1 l, D 1, 
1 Un 1° u 
(i, 7) 


Pn=0, bn Pr Ps 
1° 


af Partial Differential Equations. 


(i, 9) 


i, 
n Pi Pn (i, 
bY = 1D In—Pn ) 


(mS? +n?) 2, Pe-a Post Pn 
Lees 


n T#8 

(i, j) 
1049) 

P,=1, 4 

(i, 7) (i, 7) 


4, 
(m™ 


me? Poti Pn 


+n 
9) 
(SF Is D l, In—Pn ) 
(i, (i, 9) 


corey by 


It is evident that (a) will vanish if 


(i, 7) (i, 
Pi Pn N Pri Pn 


D4 ne? — 


The conditions (11) imply 


D4 n® >) ps Pa 


m” 


so that (b) and (ce) vanish as a result. 


Moreover, (11) also imply 


p,=0,....,0,; except that not all p’s=0; s=1,....,n), 
Hence (d) and (e) also 
Therefore, a sufficient condition that F“” may 


by virtue of (10) for successive values of the I’s. 
vanish as a result of (11). 


| 

| 

| 

| 
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replace [D™?] or {D“?} in any seminvariant without affecting its seminvari- 


ance, is that for each pair (i, j=1,....,n; i<j), 
(i, 7) (i, 9) 
™ Pr Rig 


jp + + n jp 


be seminvariants, where l, is the highest order of a derivative of [D%?] 
(or {D“?{) with respect to u, (r=1,...., 2) occurring in the original 
seminvariant. 


6. Application of the New Method to Special Cases. 


It should be noted that while we have indicated a large class of expres- 
sions which may be used in the place of [D“?] or {|D%?}, if we wish the 
simplest set of seminvariants we will set =1, n©®=0, or m®®=0, n©®=1, 

Let us now illustrate this method of computing seminvariants by the case 
of plane nets. We assume temporarily a,=b,, and take a,/A=b’, 4,/a=a’. 
Substitution of these values in (3) yields the fundamental set of seminvariants 


A, =a—2b’, B, =), C,=c +ab’+a’b—b,—b”,* 
A; =0, B, =0, C, =c’ +a’b’—b), or c’+a’'b’—a;,, + (12) 
Ay =e, By =b"—2a’, Ci’ 


which are considerably simpler than those of (4). A further advantage of 
the new process over the old becomes evident when we discuss the invariants 
of (1). Of the seminvariants (4), only two, A” and B, are invariants. Of 
the seminvariants (12), five are invariants, viz., those just mentioned and in 
addition C,=€, C}|=K (or H), Ci;=€”. That this advantage can be guaran- 
teed in general, however, is not evident. The simplicity of the new process 
might perhaps give it a better chance than the old to yield a larger number of 
seminvariants which are invariants as well. Another pseudo-canonical form, 
viz., that determined by A,/A=4a, 2,/A=3b”, yields a fundamental set of 
seminvariants which are simpler than those of (4), but not quite so simple as 
those of (12). Only two of these are invariants. 

In the case of conjugate systems of curves on a curved surface,t the new 
process (for a certain choice of 4) gives seminvariants which are very much 
simpler than those given by the regulation method. Moreover, three of the 


*It happens that C, and C,”’ may be still further simplified by use of the integrability conditions. 
+G. M. Green, AMERICAN JOURNAL OF MaTHEMaTICS, Vol. XXXVII (1915), pp. 215-246. 
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five are seminvariants, as against one out of five for the old method. In the 
case of developable surfaces,* besides yielding simpler seminvariants, the new 
method gives two invariants, while the other gives none. 

It seems evident that the method we have outlined will be of even greater 
value in the cases of higher dimensional configurations. For example, the 
simplification due to the new process for the case of curvilinear coordinates 
in n dimensions + is seen by noticing that it would replace, for every 3, 


f = 
by some single coefficient af”. 


ANN Arsor, MIcu., January 2, 1918. 


* W. W. Denton, Trans. Amer. Math, Soc., Vol. XIV (1913), pp. 175-208. 
+ G. M. Green, “The Linear Dependence of Functions,” etc. loc. cit. Cf. especially Section 7, and 
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On a Method for Determining the Non-Stationary 
State of Heat in an Ellipsoid. 


By Datta. 


Introduction. 


1. The first writer, who attempted, with some success, the problem of the 
determination of the non-stationary state of heat in an ellipsoid with three 
unequal axes, was E. Mathieu,* who showed how the problem could be reduced 
to the solution of certain ordinary linear differential equations. But he found 
these equations to be so unmanageable that he contented himself with approx- 
imating to their solutions for the special case of an ellipsoid of revolution. 
Prof. C. Niven improved upon the results of Mathieu in certain respects in an 
interesting memoir, +t entitled “On the Conduction of Heat in Ellipsoids of 
Revolution.” 

In the present paper, I propose (1) to obtain the chief results of Prof. 
Niven by an entirely different method, and (2), to show how this method can 
be applied to the case of the ellipsoid with three unequal axes, to obtain similar 
results which are believed to be new. It may be noted here mt, in Art. 6, I 
point out a mistake in Prof. Niven’s memoir. 


Preliminary Remarks and Definitions. 
2. Let the nee Nee aa of the ellipsoid be f(a, y, z) and let its 
boundary, +5 =1, be maintained at temperature zero. Then 


the required y, 2,t) is such that 


* Cours de physique mathématique, Ch. IX. } Phil. Trans., Vol. CLXXI (1880). 
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the units being so chosen that the diffusivity is 1, 


V=0 on the boundary, (2) 

V=f(2a, y, 2) when t=0. (3) 
Thus V can be expressed as a sum of terms of the form Ae~W(z, y, 2), 
where the normal function W satisfies the equation 


aw 


and vanishes on the boundary; and the constants A are so chosen that the 


initial condition (3) is satisfied, so that f(2, y, 2) =LAW. 
When a=b=c so that the ellipsoid becomes a sphere of radius a, an 


appropriate norma] function is 
S,, (ar) P™ (cos 0) cos mg, 


where 
S (x) =(— P%®(cos 6) =sin u being cos 6, 


and A is a root of the equation S,(Aaa)=0. Throughout the present paper I 
will represent by W* the normal function corresponding to the function 


S,,(Ar) P™ (cos 6) cos m@, and denote W? by W,,. 
I proceed now to obtain the functions W of various types. 


W, For Elltpsoid of Revolution. 
3. Let e denote the eccentricity of the ellipsoid, then neglecting e* and 
higher powers, the equation of the ellipsoid can be written as 


where a=a(1—4e’) and e=4e. Now assume that 
W,=S,(ar) +e=M,S, (ar) P, (cos 6), 
t=1 


M, being an unknown constant to be determined. Then, evidently, W, satisfies 
the partial differential equation (I); and to satisfy the boundary condition we 
must have 


p, + 63. M,8,(aa)P, (cos 6), 


since e’ and higher powers are neglected. This equation must hold for all 


| 
| 
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values of cos 6. Therefore, equating to zero the coefficients of the various 
zonal harmonics, we get 


Sy (Aa) =0, (1) 
M,8, (Aa) +a =0, (2) 


and all the other M’s are zero. Hence the required expression for W,, in 
terms of a and e, is 


W,=S,(ar)—4e? S,(Ar) P,(cos 6), 


S, (aa) 
where 4 is given by the equation (1). But the general solution of the equation 
(1) is known to be 

Aa=17, 
1 being any integer. Hence 
Aa=in-+t ine’. 


4. In order to obtain a closer approximation to W,, we will retain e* and 
neglect the sixth and higher powers. Thus the equation of the ellipsoid is 


where B=a(1— e’— c= + and t= Fe’. 


Let us assume that 


W.=S8,(Ar) —F S,(ar) P, (cos 6) +7 P, (cos 6), 


S,(Aa) 


N, being an unknown constant to be determined. Expanding by Taylor’s 
Theorem, we have 


and S,(ar) =S8,(aB) (oP, +7tP,) +....; 
when 
Again (P,)’= 


18 


| 
| 
| 
| | 
} 
| | 
OS, (Aa) 
| 
| 
| 
| 
| 
: = 
dS, (Aa) 
| 
i 
| 
: 
| 
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Hence we must have 


44 (18)P,. 


This equation must hold for all values of cos 6. Therefore, equating to zero 
the coefficients of the various zonal harmonics, we get 


dS, (Aa) 


08, (Aa) 
2 da 0S, (AB) 
— jes B aB (2) 

OS, (Aa) 


2 da 0S, (AB) 
S, (Aa) OB =0; (3) 


and all the other N’s are zero. On substituting the values of 6, o and ¢ in 
terms of a and e, the above equations take the forms 


Sy(a 


OS, (Aa) 

08S, (Aa) 

0S, (Aa) 


(2a) _ 2 "Oa, _ 


8, (Aa) da 


4 
4 


| 

| 

| 
| 

| 

| 

| 
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Hence the required expression for W, is 
W,=S8,(ar)— te? S,(ar)P, (cos 6) + e!N,S, (Ar) P, (cos 8) 


+ 2 eN,8,(ar) P,(cos 6), 


the values of N, and N, being given by the equations (5) and (6), and the 
value of A being determined from the equation (4). 

5. I will now proceed with the solution of the equation (4) of the pre- 
ceding article. If we neglect e altogether, the equation reduces to S,(Aa) =0, 
whose roots are given by Aa=in. Therefore let the full value of Aa be 


1,, 1, being unknown quantities which are to be determined. S,, can be expanded 
in a series containing a finite number of terms so that 


5 
n' (n’—1’) (n’—2’) (n’—3’ =) 
(n’—2') | 1 


where n’ stands for n(n+1) and (n’—?’) for (n—t)(n+t+1). Hence we 
have the following: 


2 


8, (aa) = cos int (2 — giz) cosine... 


S, (Aa) =— 35008 — cos im— im e+... 
2 
ot —2cosin+. —(1— =) cos in+...., 
| 


Therefore, substituting the above expressions for S,(Aa), etc., in the equation 
(4) of the preceding article, and equating to zero the coefficients of e* and e*, 
we obtain finally 


Aa=in+ + ine? + 135 (in?-+27)e! 
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Comparison with the Results of Niven. 


6. By neglecting e° and higher powers, Prof. C. Niven has obtained an 
expression * for the parameter A, viz.: 


im ,. 
Aa=in+ + ine®+ 105 


which differs from that obtained by me in the preceding article only in so far 
as the coefficient of e* is concerned. After carefully going through Prof. 
Niven’s calculations, I find that the mistake of Prof. Niven must be attributed 
to some inadvertence on his part, as, by repeating his process, I get the correct 
result. 

From the identity of the values of A, it follows at once that my expression 
for W, and Prof. Niven’s must be identical. For, as is well-known, for the 
same value of 4 there can not be two different solutions of the equation 

ow, FW, 


both of which vanish on the boundary. 


+v’W,=0, 


W,, For Ellipsoid of Revolution, n>0. 


7. Let e* and higher powers be neglected, so that the equation of the 
ellipsoid is the same as in Art. 3, viz.: r=a}1+eP,(cos@)}. Then assume 
that 

W,=S, (Ar) P,, (cos 6) +eL’H,S, (ar) P, (cos 6), 
#=0 


where the H’s are unknown constants, to be determined, and >’ refers to all 
t=0 


the values of ¢ except i=m. Thus it is evident that W,, satisfies the partial 
differential equation (I), and it remains to find the values of the constants H’s 
so as to satisfy the boundary condition. Since «? and higher powers of e are 
to be neglected, we have, on putting r=a(1+éP,), 


8, (ar) =8, (Aa) p, 


also, we have 


where 
(ntl) (m+2) _ a(n) 
(2n+3) (Qn+1)’ (2n-+3) (2n—1)’ 
n(n—1) 


(2n+1) (2n—1) 


* Loe. cit., p. 145. 
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Hence we must have 
OS, (Aa) 
0a 


This must be true for all values of cos 6. Therefore equating to zero the 
coefficients of the various zonal harmonics, we get 


0=S,(Aa)P,+ea Py. 


0S, (Aa) 
+ a 2°82) (2) 


and all the other H’s are zero. Thus the unknown constants are determined 
and the required expression for W,,, in terms of a and e, is 


dS, (Aa) 
da 
Si42(Aa) 
OS, (Aa) 
a 
da 
S,-2 (Aa) 
A being a root of the equation (1). 


a 


W,=S,(ar)P,—+e S42(Ar) Pi 


But, in terms of a and e, this equation is 


n?+n—1 2a OS,(Aa) 
(2n+3) (2n—1) 


(aa) — 0. (4) 


Therefore, if x be a root of the equation S,,(z) =0, the corresponding solution 
of (4) is given by 


n?+n—l1 


(2n-+8) (2n—1) 


Ws For Ellipsoid of Revolution. 
8. Let e* and higher powers be neglected, and assume that 
(Ar) P™ (cos 6) cos mo+e (ar) P™ (cos 6) cos 
t=m 
where 2’ refers to all values of ¢ from m up to infinity except the value t=n, 


t=m 


and 7” is an unknown constant to be determined. Then it is evident that W* 


| 
| 139 
| | 
j 
| | 
i 
i 
{ 
| 
H : 


140 Datta: On a Method for Determining the 


satisfies the partial differential equation (I). Now, putting r=a(1+eP,), we 
get 


S, (ar) =S, (Aa) +ea 


where 


and 


(2n+3) (2n+1) (2n+3) (2n—1) ’ 


(n+m) (n-+-m—1) 
(On+1)(2n—1) 


Hence, from the boundary condition, we must have 


X cos I? S,(Aa)P™ cos md. 
t=m 


Therefore, equating to zero the coefficients of the various surface harmonics, 
we get 


8, (2a) —o, (1) 
0S, (a 
+ a4) (2) 
S,, 
a4) pp (3) 


and all the other J?’s are zero. 
Thus the required expression for W™ in terms of a and e is 


dS, (Aa) 
W2= |8,(ar) (cos 6) — 0) Bs 
L Sn42(Aa ) 
(Aa) 
Oa 


where 4 is a root of the equation (1). But, expressed in terms of a and e, 
this equation becomes 


8,(aa) — (n?-+-n—1) +m? (Aa) (4) 


(2n+3) (2n—1) da 


Therefore, we obtain 
Aa=x {1 


(n?-+-n—1) +m? 
(2n+3) (2n—1) 


corresponding to the root x of S,(#)=0. 


| 

| 
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W, For Ellipsoid with Three Unequal Axes. 


9. Let e, and e, be the eccentricities of the two principal diametral sec- 
2 2 2 
tions of the ellipsoid =1, a>b>c, by planes passing through 
the major axis. Then, neglecting the fourth and higher powers of e, and ¢,, 
- the equation to the ellipsoid is 


(P»—P.) + (e2—e?) Picos 29], 
4. cos 29), 
where y=aj1—<(ei+e)}, a=(eit+e), and 
Now assume that 


W,=S,(ar) + R,, (cos 6) cos me, 


t=m m=0 


where #,, , is an unknown constant to be determined and >’ refers to all values 


t=m 
of ¢ from m up to o, except t=0. Then it is evident that W, satisfies the 
partial differential equation (I). To satisfy the boundary condition, we must 
have 


0=S, (Ay) (e,P.+6P3 cos 2) +2 Ry PP (cos 6) cos m@, 


the unknown constants R’s being assumed to be of the same order as ¢, and e&. 
Therefore, equating to zero the coefficients of the various surface harmonics, 
we get 


S,(ay) =0, (1) 
Ry» (Ay) + ey =0; (3) 


and all the other R’s are zero. Thus the R’s are determined, and we get 
finally the required expression for W, to be 


(Aa) 


a 
+ 4 Ps (cos 6) cos 29, 


| 
| 
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| 
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where A is a root of the equation (1). But the roots of the equation (1) are 
given by 

Ay 

Hence, in terms of a and the eccentricities, we have 


Aa=in-+ + 3) in. 


Conclusion. 


10. I conclude this paper by pointing out that the results or Arts. 3-9 
admit of extension and generalization in various directions. For example, a 
procedure similar to that of Art. 9 will give us W® for the ellipsoid with three 
unequal axes. Also, denoting by W™, the normal function corresponding to 


S,,(Ar) (cos 6) sin mq, 


it is obvious that, for the ellipsoid of revolution as well as for the ellipsoid 
with three unequal axes, W™ can be obtained in the same way as W”. 


UNIVERSITY OF CALCUTTA, 1917. 
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Nilpotent Algebras* Generated by Two Units, i and j, Such 
That i’ Is Not an Independent Unit. 


By Guy Watson 


I. Introduction. 


The problem of referring all hypercomplex number systems to a relatively 
small number of typical forms was suggested by Hamilton,t but with the 
exception of DeMorgan’s discussion of double and triple algebras, nothing 
much was done till Benjamin Peirce} worked out all algebras of deficiency 
zero and one. Starkweather$§ worked out algebras of deficiency two. He 
showed also that algebras of nm units could be obtained from those of n—1 
units. Cartan || using the characteristic equation developed the semi-simple 
and the nilpotent sub-algebras, and showed the possibility of representing 
every algebra by means of units with double character. Taber { reestablished 
the results of Peirce and extended them to any domain of rationality for the 
coordinates. Wedderburn** and Vogherat+t made an advance in the treat- 
ment of the hypercomplex algebra by basing their work on the conception of 
invariant classes of numbers in the algebra. 

Besides this direct line of development there have been two others. The 
first is by means of the continuous group, the second is by using the matrix 


* This paper considers only linear associative algebras whose coordinates are taken from the field 
of ordinary complex numbers. 
7 “ Lectures on Quaternions,” Preface, pp, 29-31. 
+“ Linear Associative Algebra,” AMERICAN JOURNAL OF MATHEMATICS, Vol. IV (1881) pp. 97-192. 
§ “ Non-Quaternion Number Systems Containing no Skew: Units,” AMERICAN JOURNAL OF MATHEMA- 
Tics, Vol. XXI (1899), pp. 369-386, 
|| «Les groupes bilineaires et les systemes de nombres complexes, Annales de la Facultedes Sciences 
de Toulouse, Vol. XII (1898), B. pp. 1-99. 
q ‘On Hypercomplex Number Systems,” Am. Math. Soc. Trans., Vol. V (1904), pp. 509-548. 
** «On Hypercomplex Numbers,” London Math. Soc. Proceedings, Series 2, Vol. VI (1907), pp. 
77-118. 
tt “ Zusammenstellung der Irreduziblen Komplexen Zahlensysteme in sechs Einheiten,” Denkschriften 
der Math. Nat. Klasse der Kais. Akad. der Wiss. Wien., Vol. LXXXIV (1908), pp. 269-329. 
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theory. The first method was used by Scheffers,* Molien,t and Study,{ the 
second by Shaw § and Frobenius. || 

Shaw regards all associative numbers as belonging to an associative 
algebra of an infinite number of units, the “associative units,” 4,,,, which are 
elementary matrices. Each associative number is a linear combination of 
these units, so the theory of linear associative algebra is the theory of these 
associative units. He shows that the presence of a modulus is not necessary, 
thus making the methods particularly applicable to nilpotent systems. He 
proves that the equation of an algebra determines all the units but those which 
form a nilpotent system, and consequently to get all linear associative algebras 
we must first determine all nilpotent algebras. Benjamin Peirce] was the 
first to recognize the importance of nilpotent algebras. Furthermore, algebras 
of order » may be found without first knowing those of order n—1. By 
selecting a base and adjoining to this a nilpotent unit, an ever-increasing 
system of nilpotent algebras may be determined. 

The simplest case of the Shaw canonical form ** of a nilpotent algebra is 
that in which there is only a single generator, 7, whose u,-th power vanishes, 


Three lemmas concerning polynomials in this nilpotent have been introduced 
for use in the handling of the next simplest case, namely, that in which there 
are two generators, 7 and another nilpotent unit, 1, whose square is not an 
independent unit in the canonical form. The unit j is such that j“~' 0, 7“=0 
and ij""=£0, ij#=0, where uw, and uw, are multiplicities of j relative to j and 
to respectively, and The expressions i, where 0<s<y,, and 
0<t<y,, are the independent units of the system. The algebra is therefore 
of order u,+“.—1, or if we insert a modulus y, of order uw,+u,. The u, power 
of every number must vanish, hence the deficiency in this case is u.. There is 
at least one hypernumber which does not vanish for a power lower than 4, 


* “ Zuriickfiihrung complexer Zahlensysteme auf typische Formen,” Math. Annalen, Vol. XXXIX 
(1891), pp. 293-390. 
¢ “ Ueber Systeme héherer complexen Zahlen,” Math. Annalen, Vol. XLI (1893) pp. 83-156. 
£ “Ueber Systeme von complexen Zahlen,” Gétt. Nach. (1889), pp. 237-268. ‘ Complexe Zahlen und 
Transformationsgruppen,” Leipzig Berichte, Vol. XLI (1889), pp. 177-228. 
§ “Theory of Linear Associative Algebra,” Am. Math. Soc. Trans., Vol. IV (1903) pp. 251-287. 
“On Nilpotent Algebras,” Am. Math. Soc. Trans., Vol. IV (1903), pp. 405-422. 
|| ‘ Theorie der hypercomplexen Grossen,” Berlin Berichte (1903), pp. 504-537, 634-645. 
q Loc. cit., p. 118. 
** Loc. cit., p. 406. 
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and there may be several such. The base unit i is chosen such that ij“=0, 
ij“—! 0, the various products 7, 17*,...., do not contain terms in 7 alone, 
such as aj”, and yu, is such that there is no hypernumber 7’ satisfying the con- 
ditions on for which 7’j%=0, $0, In other words, is not the 
product of any hypernumber into j, and we do not have for any power of J, 
ij°=bj'. These characters of «1 and j are essential. The two sequences 
jy and 4, 17, ...., constitute the two shear regions of 7. 
In the cases Peirce considered, that of deficiency zero contains no 1; that of 
deficiency unity must have u.=1, hence ij=0. In Starkweather’s types 
where only two generators enter, u4,.=2, hence 7 is a unit but 77?=0. In the 
types herein considered the deficiency does not play any role at all. The 
investigation is along the line of Shaw’s construction by generators and not 
by classification by numerical invariants, other than those entering the equa- 
tions of condition. The associative units are used, though not indispensable, 
for convenience in making reductions. 


II. Lemmas Concerning Polynomials of a Single Nilpotent 7. 


Lemma l, If j is a nilpotent number (j“=0), and if F(j) and G(j) are 
polynomials in j, a quotient Q(j) can always be obtained Gen provided F(j) 
does not contain a term of lower order than does G(j). 


Proor: The lemma will be proved if we can find a polynomial, Q(j) of not 
more than yu, terms such that F(j) =Q(j)- G(7), or 


ther? = (Got +927? +--+) (Got 


Equating coefficients of corresponding powers of 7 on both sides since they 
belong to independent units, 


These equations can not be solved if F(j) starts with a lower power of j than 
G(7), that is, if 


+0, 
while If, however, 
while FO, O<k<h+1, then the first h—k+1 
of the above equations are identically satisfied, the next k equations give 
and from the remaining u,—h—1 equations, q,, 


Gu, arbitrary. 
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Lemma II. Consider the equation L(j)-M(j)=0 where 
L(j) M(j) = My + + M99? + +m, 

(a). If 1,460, M(j) vanishes entirely. 

Then 

omy =0 .°. m=0, Lom, +lm=0.°.m,=0, 

vanishing terms =0 

(b). If 4==1=....=1,=0, 14,40, then 
M (9) (J) =0, where t+ +i, and 0. 
Applying (a) to the equation M(7)j"*'- Z,,1(7) =0 gives 


where M,(j) =my+myj+mypj?+....+m,j"+ any terms whatever up to j“—. 
Lemna III. If C(j) is a polynomial in j with non-vanishing constant 
term, then j can be expressed as a polynomial in (jC). 
Proor: Let C(j) #0, then 


(jC)? + 2a,a,7° + (jC)? 
And since a), a2, a3,...., a! 0 by hypothesis the determinant of coefficients 
does not vanish, and we may solve these equations for j, 7°, 7°*,...., j#7? in 

terms of powers of (jC). 


III. Eapressions for the Generators. 

From a nilpotent system choose any nilpotent expression 7 to be a unit, 
called the adjunct unit, such that 7 is a number which has as high a non-van- 
ishing power as any number of the system, and a set of expressions (in this 
paper one, say 7), called the base. Then it is known that any hypernumber of 
the system is linearly expressible in terms of the u,+y,.—1 independent units 


and in this order the product of any two units 77" and 77” is linearly expressible 
in terms of the units which follow* the unit ij**?—. 


*Shaw, Trans. Am. Math. Soc., Vol. IV (1903), pp. 406-410. 
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It has been shown + that any number is expressible linearly in terms of 
elementary matrices, A,,,, called “associative units,” + where 
dSw=lifs=r, Vifs#r, 
c=1 if t+?’ satisfies the above condition, 0 if ¢+¢’ does not satisfy the condi- 
tion, wu, and uw, are multiplicity numbers for r and s respectively, determined by 
the power of 7 which equals 0, and the power of j in ij” which causes it to 
vanish, that is here j7“=0, 0 and ij”=0, 17% * $0. 
The frame of the nilpotent system formed upon the base 7 and the adjunct 
unit j is:* 
We may indicate these expressions more compactly by introducing certain 
symbolic polynomials which will take care of the third subscripts. Thus we 


write 


where A .... + and is a nilpotent such that 6*=0, 
Aco9=Aco;41. In a similar manner we may write 


where + +++ and 6 is a nilpotent such that 
5“=0, AywS=Awi41. We need to notice then that we may also indicate in the 
Same way a series of terms in (,,, 


= yy Agu + 


where + 


Hence in this notation, if we omit A and write only subscripts, 
i= 210+ 22044126" 

where 

If 4:=4., 6,>=b,=0 for 121 can not appear in i. We also use the same 

symbols for polynomials in 7, the context always showing the use. 


* Shaw, Trans. Am. Math. Soe., Vol. IV (1903), pp. 251-287. 
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The general form of the product 72 is necessarily 
+ in which e,=e,=0. 
Tueorem I. If C1 by a proper choice of 1 the product ji becomes 
ji=tjC(j). 
Proor: Let X=xj+2,j?+...., then 
j(i—X) =jC+H—jX= 
If we can determine X so that E+7X(C—1) =0, we may set 7’=i—X, whence 


jv’ =v'jC, which is the form desired. We must have then 


g(1—C) 1—C 

By Lemma I, X is determinable provided the numerator does not start 
with a lower power of 7 than the denominator. This is evidently the case 
whenever 1—C starts with a term free from j, but when C=1+ )7C,, where 
C, starts with a term free from 7, then the numerator must start with a higher 
term than j*~*. To show that in this case the desired situation always exists 
we proceed as follows: 

Let us abbreviate the writing by setting temporarily 

V=YA+B, pi=yC+j"E,, 
where EF, starts with a term free from 7. Then we can easily see that if P(7) 


is any polynomial in J, 


P(j)i=iP(jC) when C1, 


where the division is purely formal and always possible. 
We have now j:?=j7ijA+jB=i77AC+ terms in j. But we have by 
changing the association, 


C(jC)—1 ... + terms in j 
C(jC)—1 
C(jC) +ij"C ii (50) | + 


Since the terms in ij? must be equal, we have 


~ 
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Now suppose that C=1+ °C, where C, starts with a term free from j, then we 
have, since we may drop the factors that are invertible, by Lemma I. 


Ci 


Now let A=j”A, where A, starts with a term free from j, then the lowest 
terms from the two expressions that are multiplied by 7 are 


+I" 


where a, is the first term of A,, c, of C,, and e, of H#,, none of which are 
zero. Hither these two cancel each other, which can happen only if 


a+y+2=w, and hence x<w—2; 
or else each is not lower than u,. But we know that x<u,—1, and hence if 
we have 


In any case then we must have x<w—2. Therefore whenever C1, 


X(j) is determinable so that jidgGy. 


Moreover, we have X=j"~"(X’), and if /#=i—X, 
Hence 7’ is of the same character as i, and the transformation is permissible. 
In other words, whenever C1, we may take E=0 and write, 


j=111+4 2200. 


We proceed to investigate first the case C=1. 


IV. Case I. C=1. 


If C=1 the generators have the form | 


Since, however, in this case 
| 


we may set i’=i—4jA, whence 

v? =?—[ijA-+ terms in powers of j] = terms in powers of j. 
Consequently, we may treat the case as if Ad=0, until we derive the types. We 


then return to the equations in i, 7. The generators now become 


t= 210+128"-“B, 
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From these we form the following products: 
47 = 2114 220% 41254 = 21 n+ 
7G (j) =216G (6) +220" “FE BG (6), 
ji= 2114126" 
(i Co + — +1) e,5+....], 
The last two terms must vanish, so we have 
—0, (1) 
P= 226% 
This gives rise to the equation 
(ty — Le) Dot (tr — +1) ... .] =0. (3) 
The solutions of these three equations are in three classes: 


EO. 
3. 
For Class 1 equation (1) is satisfied. Equation (2) becomes . 
This equation is satisfied unless a=0 and then H=dE,. Equation (3) becomes 
(u.+a)bo+....]=0, which is satisfied whenever (2) is. 
For Class 2 let where e,>-0. Then (1) becomes 
Equation (2) becomes €,....]=0, which is satisfied 
unless v=2u,—y,, in which case E,=dH,; and H=6""'E,, therefore v >2u.—. 
Equation (3) is satisfied with (2). 
For Class 3 all three equations are satisfied. We have then for this case 
the following types: 
12. m=24,, H=dE,, 
3. E=0. 
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These types are now written out as follows: 
11. a>0,H The generators for this type are 


4=210+4+126"t*B, 


From these we get the following products: 

47 72 =1124+ 2224125444127, 

13%. =2ue, E=6EF,, E, #0: 

4 =210+126"B, 17=2114+125"*'B, 

2114-126" B 41108, (7), 

=2104126""B, , 

47 = 2114120 7? F, , 

ji = 2114-12581 B +1154 = + (9), 

—2294-4B + 116“-“B (7) 
3. E=0: 

=2104126"-“B, j7=1114221, 17 

ji = 2114-12041 B=1j, 7°=11n+22n, 

?2—1164-“B + B— B ( 


Type v= (’)?= 
1 
3 j@**B(j), 
or in terms of 1 and /, instead of 7’ and j. 
Type 
1 if, (j4)'’—47°A? +9" (4), 
2 (9), (7) (9), 


when (jA)’ is the formal derivative as to 7 of the polynomial jA, and A is 
arbitrary. In case A starts with a term of order u,—y,—1 or higher, «’ may 
be used for 1. 
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V. Case II. A=0, C#1, E=0. 
If A vanishes we have the defining equations: 
4 =210+120"-"B, j=111+226C. 
Then and iP(j) =21P (6) +126"-“BP (dC). 
Again ji (j) —126" “BC - C (6C) +1264“ 
(1) 
, whence P(j) =11P+22P(6C), 
“B, 
[B—C*-”B (6C) =0. (2) 
Let C be written C=c,+6'C, where c, may be 0, but c,~0 and 0<t<yu,—1. 
Substituting in equation (1) 
— +6°C,) ) ] =0, 
or +1B (0C) ] =0. 
By applying Lemma II we have either 
11. and 
or 13. 6=l ead 
Now substituting in equation (2) the results just found, we have for 11: 
s[1— =0. 
Therefore by Lemma II, 
111, b,,1.=0, and B=0; or 112, c#=1, b,, 
For the case 12, B,—C"—*1B, (6C) ] =0. 
Hence either 121, for c.=1, 6“~by)(u,—t—1)C,=0, for which we have either 


1211, m>m, 
or 1212, 
or 122, for (—1+ 6'C,)"—*"B, (6C) ] =0, 
for which we have 

1221, u,—1 is even, B,=6'B,, t.e., B=0""D,,_1, 
or 1222, B,=0""B,, B=0, 
or 1223, w,—t—1>., By #0. 


q 
| 
| 
| 
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Consequently when A=0 the following seven types arise: 
A=0, C=c,+6C,, 
B=0. 
112. 
1211. q=1, a> m- 
1212. aq=1, 
1221. q=—1, u,—1lis even, 4. 
1222. c=——1, B=0. 
1223. Co= —l, wm, B, +0. 
We now work out these seven types obtaining the expressions for the 
products ji and 7. 


111. C=c,+6'C,, @&-1, B=0. For this the generators have the form 
4=210, 7 =111+220C. 
Then ij=211, and ?=0. 
3. 


Now i j=111+220C, ij=211, ji=210C=ijC(j), 


1311. q=i, 
§ 111422142260, ij" = , 
ji = 211 +. 21010, = ij +. 
j*=11n + 226°C" =11n + 226"(1+6'C,)", 

1212. q=1, B=d""B,, 
j=111+2214+ 226°C, =111+4 220, 
+220"'B, = j"'B, (j). 

1221. is even, c=—1. 
=210+128"—"b, =1114221(—1) +226°C,, 17=211, 
ji=211(—1) +2100, = 
118") + 22641 = 

1222. q«=—1, B=0. 
§=210, j=111+221(—1) +220C,, ij=211, 
ji= 211 (—1) +2100, = —ij  #=0. 
20 


+ 
i 
+ 
| 
i 
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1223. B=d“"*"B,, B, +0. 
7 =111+4221(—1) +226+C,, 
47" 
jt=211(—1) + 216°C, 4-126" B, = —ij (7), 
(4). 
Type 
1 49 0, #1, 
3 (9), ge (9), > lle» 
6 (7); 0, 
7 4—t—124,, B, #0. 


VI. Case lI. C#1, E=0. 


In this case, which is the most general, we first form the following 
products from the generators: 


j=111+220C, 

r<u.—l, 

2101412547 | 

iP (j) = 210P (6) + 2264 P(6C) + 126"-"BP(dC), where P(j) is a 


polynomial in j. 


P(j) =11P(5) +22P(6C), 
ji=210C 

=ijC (j) - C(0C) —126"-“*'BCC (0C), 
P(j) (6C) +2204 P(6C) +125"-"BP(5) P(jC), 

A(j) + + 2202424 4. B 

—22PACA (0C) (6C) —126"-"BCA (0C). 


Since the terms in 22() and 12() must vanish in both the products ji and 7, 
we have four equations: 


Before reducing these we need a second theorem which we proceed to deduce. 


(1) @AC[1—C(6C)]=0, 
(2) —C-C(8C)]=0, 
(3) A(8C)]=0, 
(4) @A[A—C- | 


‘i 
| 

| J 
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If A(j)=j°A,(j), OSs<m—1, then “B(j), where 
A440. By Lemma III we may now express a polynomial in j in terms of 
one in jC if C=c,+jC,. Let A,(j)=A;(jC) and determine 7’ from i=7' A; (j) 
where a, 0, then 

But using the formula P(j)i=iP (jC), 
Aj(j) +4 (jC) =14, (9), 
Aj(j) (7) Ay (7) (9). 
Thence (9) Ai (9) +9" “B (9). 

Tueorem II. If where c, 0 then j=P(jC) and A;(j) 
is known, and invertible. Hence we may choose i so that the term 1j**'A;(7) 
becomes simply 17°*', that is, we may take A,=1. 

As a consequence of Theorem ITI we have two sub-cases, 

A(j)=7', 
t>0, Co #9. 
Case ITI,. The generators may now be written 
i= 210+ 226414126" "B, j=111+220C, 
and we shall have 
ji=yC(j), 
where the following four equations, which are the reduced forms of equations 
(A), must be simultaneously satisfied : 
(1)’ #&?C[1—C(6C)]=0, 
(2)’ - C (dC) ]=0, (A)’ 
(3)’ 
(4)’ + [| B—C*-“B (6C)] =0. | 


Consider equation (1)’, 
[1—C(6C) ] =0. 


This leads to two classes 3 and 4. 


3 Co FO. 
4 s+2<u,, then since or C(6C) 
= 1+ Cag FO. Co 
or Cot 6 [Cy (OC) ] +6'C,] 


| 
| 
| 
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Therefore and t=u,—s—2+r since and C(6C)=1+6C,, 
t<,. Consider next equation (3)’, 
For class 3, C=O+6'C,, 0, #0 this becomes 
6“—16, [1—C”-1] =0, 
from which we have either 
31 b)=0, B=6B,, i. 
or 32 


For class 4, s+2<u,, t=yu,—s—2+r equation (3)’ becomes 
This is satisfied unless r=0, then B=6B,. Consider now equation (2)’, 
C(dC) ] =0. 
For 31 s+2=,, ce 1, C=o4+8C,, Put B=sB,, 
Then .C (dC) ] - C(dC) J. 
311 If then B=d"~"d,,,, or B=0 if b,,,=0. 
312 If ci=1, that is c=—1 and is odd, then 
.°. v2=m—t—1. 
For 32 ce——1, then 6“-“t1B[1—C - C (dC) ]=0, 
BT 1 —(¢)+6'C;) (C9 +8'C'C, (SC) ] =0. 
321 If 
322 If c=1, B=d"*"B,, and if ¢=—1, is even, that is p, 
is odd. 

For Class 4, s+2<., equation (2)’ becomes 
Bl1—(1+6'C,) (1+8'C'C,(6C) ] =0, 

We consider finally equation (4)’, 
+ [ B—C*-”B (6C) ] =0. 

For 311 B=d”-'b where b may be 0, s+2=m, ch? FOF cy, 
then —'b[1— (cq) + =0, which is satisfied unless u,=p, since *=0. 
Therefore we have also 

3111 t.e., v>u,—1, and 

3112 w=", ch’~'=1, but this is impossible since it is a contradiction 
to 311. 


* We write 6 in place of by,—1 in what follows, 


| 
| j 
| 
| 
| 
| 
| 
| | 
] 
| 
| 
| 
| 
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For 312 c=—1, is odd, B=6°B,, by v>u,—t—1, equation 


(4)’ becomes 
(—1+0C,)"-"*"B, ] =0. 


We have then either 


3121, is odd, and Me, 1.€., V2 OF 
3122, is even. 


Now [B.—B,(6C) eee ] 
becomes (1—C)+....]=0, 
and +1 7 fe, U2 


For 321 B=d""b,,_,, 1=1, equation (4)’ becomes 
[1— +6'C,)"*] =0 


which is satisfied since #*=0. If If w=, the [] causes 
the expression to vanish. 

For 322 ci=1, B=d"*"Bs, is even if c=—1 and (4)’ becomes 
Bs — (co + (6C) ] =0, for which we have 


3221, —....]=0, 
.*. 6*—'B;(0C) =0, which is satisfied unless Then 
32211, ps. 
32212, Bs(0C) =O6B,, t.e.. B=O*“B. 
3222, c=—1, u.—1 is even, then 


B,— (—1) "1B, — (0C) +. =0. 


32221. If u,—t—1 is even, 6“ +....]=0. 
.'. Bs=6B,, which makes B=6""D,, 
32222. If u,—i—1 is odd, B,=6*t'B,, which makes B=0. 
For Class 4 s+2<m., q=1, t=u.—s—2+r, 
Equation (4)’ becomes 


(1 —C**!] [| (6C) ] =0, 
(14 (80) ]=0, 
+.... +0°C,B,(0C) +....] 
.*. 6B, (0C) =0, for which we have either 
41, OF 


0, 
0 


? 


} 
+ 
| 
| 
| 
| 
| } 
| | 
i 
| 
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We have therefore in Case III, the following eleven types: 


a1. a>m, Ba 
3111. B=, $#1, #1. 
3121. ch 11, is odd, B=0’B,, is odd, 
V>Me—t—1, 
3122. q=—1, w—1 is odd, B=0°B,, is even, 


+ v2 


32221. c=—1, u,—1 is even, u,—t—1 is even, 
(32222. cy=—1, u.—1 is even, u,—t—1 is odd, B=0. 


C=1, t=m—s—2+r, B=6t'"B,. 


41. 
42. m=, G%=1, t=u4—s—2+r, 


These particular types are written out as follows: 


311. 


Then 


3111. 


Then 
3121. 


3122. 


7=111+220C. 
s+2=m, B=0, G1, 
4  7=111+220C. 
G=—1, 1, is odd, B=6’B,, is odd, 
U>We—t—1, 2 
4 , 
j =111-422(—1) +2260, =111-+4+2200, ij>=21n +125" 
ji=211 (—1) 1B, 4 = — ij 
= 210" 4 +B, +1196" "B, 
+9" (9). 
"$1, w—1 is odd, B=OB,, is even, 
V2 v>u.—t—1. 
4 , 
j =111+4 221(—1)+ 226°C, =111+4220C, ijn 
jt=211(—1) +216°°C, = (7), 
21014-1186" B, B +B, (7). 


| 

| 

| 

| 
| 

| | 

| 

| 
‘ 
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32211. 


32212. 


32221, 
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4 =210+ j=111+220C, i17=211, 
G=1, m>m, 
j =1114+2214+226'HC,, ij? 
ji=211 +128" B, 4.2100, +ij*1C,(j), =11n +2261 +0'C,)", 
210714 4 119" 4 

+ 220" = ij" * + (9). 
m=m, #+2=—m. 
211 +2100, B, = ij +710, (4), 
P= 226" + 210 
+ (j). 
Cy=—1, is even, w,—t—1 is even, B=6""D, s+2=p. 
ji=211 (—1) 21010, = =11n + 226"(—1+6°C,)*, 
= 216+. 226-16 


Since 6“~'=0, unless u,=,, these last two terms vanish separately. If 
they destroy each other since 


32222. 


41, 


42. 


Co=—1, is even, u,—t—1 is odd, B=0, 
=210+220"", j7=111+221(—1) +226°C,, ig=211, 
211(—1) +210 C, = — P= 210 
@ 7=1114 2214 220C. 
= + 226 
jt =210C + 226°C +126" (7) (7). 
+1164" B, + 126-44 B + 
@ j=111+220C, 
ij" = 21 n+ 4 B,C", 
jit =216C + 226°C (7) =ij (J), 
= 2161+ 22647 4 + 116°t?-"B, 
= 


| 
| 
| 
| 
| 
| 
| 
| 
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Type 


2 "Bs (9), 
8 i", is even, is odd, 
9 we is odd, is odd, 


Case III,. ¢=0. 
C=6'C,(0), 0<t<u.—2. We go back to equations (A) which now 
become, 
(1) 6*?4[1—C(0C)]=0, for which we have either 
1. t=u,—2, C=6~*c, where c is written for c,,., or 
2. 
(2) where b may be zero. 
(3) d&-"1B[A—C-A(0C)]=0. This is satisfied by the result of (2). 
We now use the results of equations (1) and (2) in (4). For 1 t=u4,—2, 
B=d""b, equation (4) becomes 
(4) (0C) =0, 
The three possible solutions of this are 
11, A=6""A, where m=p,/2 if uw, is even, 


m="—" it uy is odd, 


12, A=07 A;, B=0, uw is even. 
13, A,, ah+b=0, is odd. 


For 2 C=6'C,, B=d""b, A=0***4,, (4) becomes 


(4) (6C) ] +0"—1b=0, 


This reduces to 


| 
| 
| 
| 
| 
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because 4. €., or 
f2a—t—220, t<u.—2, which is true for this case. 
Now the possible solutions of 
424 
are 21, “4,>:, which gives rise to 
211, 


and 212, 2t>u.—2, 4,;=6'-"4,, m= if u, is even, 


2 


and 22, w4=“.=u, which gives rise to the following: 
221, 2¢<p,—2, b=0, 
222, 2t>u,—2, b=0, is even, 
223, 2t>u—2, As; a, +b=0, is odd. 
Therefore we have the following eight types: 
11. 4=6""4,, B= C= m= if is even. 


m= if u, is odd, 


m= if is odd. 


12. 18 even, A=0'F A,, B=0, C=*e. 

13. is odd, A,, a2, +b=0. 
211. tHe, —2, B=, COC, 

212. > le, 2t < C=0@C,, 


5 if uw, is odd. 


m= fs if uw, is even, m= 
221. A=**4,, B=0, C=0C,. 
222. is even, 2t>u—2, A=02*4,, B=0, C=0'C,. 
223. is odd, 2t>u—2, A=OT C=O'C,, a +b=0. 
These eight written out give the following: 


m= if is odd. 


=111+4+226c, 
Hence ij=211, #=210"4, +118" b=ij"A,(j) 
12. is even, A,, B=0, 

=210+2207A,, j=111+226"c, ij=211, 


| 
| 
| 
| 
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13, is odd, 4,, B=d"—'b, a+ b=0. 
=2104+220° 7=1114226"c, ij =211, 
210" 
=216°F 422054118" b= ij A,(j) + 
ji= 210°C, = (j), 210", 4-118" = 
m= if is even, m= if u, is odd. 
221. 2¢<u,—2, B=0, C=FC,. 
=2104+22614,, j=1114+220'HC,, ij=211, 
ji=210'10, P= 21014, = 
222. wis even, 2t>u—2, B=0, C=0'C,. 
=2104+2207A,, j=111+220C,, ij=211, 
223. wis odd, 2t>u—2, Al, B=S""b, C=O'C,, a2+b=0. 
§ =2104+220°F j=111+220C,, ij=211, 
ji= 2160, (7), 
Al 4-118") +. = ij "FA 
Type ji= f= 
2 ij? As(4), is even, 
ijt A,(j)+bj*, wis odd, b+a3,=0, 
if is odd. 
(7), (J); 2tSu—2, 


2t>u—2, u is even, 


F423), 


| 
if 
if 
| | 
| 
| 
| 
| 
7 
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Class 
I 


II ; 


III ; 


IV 


VI ; 
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VII. Complete List of the Types Arranged According to the Form of the 


Product ji, with Sufficient Conditions to Distinguish Each. 


Type 
1 ¥, 
(9); 
3 (9), 
4 
6 
7 
10 
12 (7), 
13 cig 
(7), 
15 (9), 
16 (9), 
17 
18 (9), 
19 (9), 
20 —ij+1j'"C, (9), 
21 (9), 
22 
23 
24 
25 
26 
27) 
28 
29 
(9), 


it (A), 
(7) ’ 
(9), 
""B, (9), 


ij 
ijt (7), 
ij? 

0, 

+ 

0, 

(7), 
ij™Ao(j) +5", 
ijt As(J), 

A,(j) 
(7) 
ij" A, (j) 


ijt Ae(j), 
AL(j) 


Conditions 
E=0. 
$0, 
ta, a>0, E+0, 
[41> Le, 
[41> 
[41> 
#1, 
OF le, 
CF, 
ati, 


b 


is even, is odd, 
B,(j) #0, 

is odd, is odd, 
is odd, 

= 18 even. 
> le 

is even, 

b+ai=0, uw is odd, 


m= if is even, 
m= if is odd, 
2t <u—2, 


2t>u—2, w is even, 
2¢>u—2, is odd, b+ag=0. 
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These types when arranged according to the form of the product ji fall 
into six classes. 

Class I, made up of type 1, is the only commutative type of the entire set. 

Class II, types 2-4, have for the product 71 the term 77 and powers of 7, 

Class III, types 5-10, have ji=1j-+77't'C,(j) where ¢>0 and C,+0. 

Class IV, types 11-15, differ from those of Class III in that 

where ci 1. 

Class V, types 16-22, are essentially different from those of Classes III 
and IV since ji=—1j+17't'C, (7). 

Class VI, types 23-30, have the product ji starting with at least 1)”. 

We need now to see that the types of each class are essentially distinct. 
The conditions given, having arisen in the determination of these types, are 
sufficient to show this. For example, we might think that type 18 and type 20 
would be the same if in type 20 b=0. Although this would make the products 
ji and @ have the same expressions for the two types, they would still be dis- 
tinct, for in 18 u,—? is odd and in 20 it is even. Types 13 and 15 have the 
same form for ji and 7, but they are essentially distinct because in 13 w,>,2, 
which is not necessarily the case in 15, and furthermore in 15 c}-'=1, but the 


only limitation on c¢ in 13 is c}61. Similarly it can be shown that the thirty 
types are all distinct. 
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